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Data  were  obtained  on  Clinton  P.  Andefson  Los  Alamos  Meson 
Physics  Faciiit/  Energetic  Pion  Channel  and  Spectrometer  for  elastic  and  in¬ 
elastic  pion  scattering  from  ground  state  nuclei.  States  observed  here  in¬ 
cluded  those  of  0.0,  5.27, 6.32, 7.16, 7.30,  7.57,  8.31,  8.57,  9.15,  9.76,  9.9, 10.7, 
11.3,  11.9,  12.5,  12.9,  13.1,  14.1,  14.4,  14.6,  15.0,  16.5,  16.9,  17.2,  17.6,  18.3, 
18.7,  and  18.9  MeV  excitation  energies.  Angular  distributions  were  obtained 
for  scattering  at  angles  from  25^^  90^in  5'^increments  with  an  incident  pion 
energy  of  164  MeV.  Optical  model  anal^es  of  the  elastic  (0  MeV)  angular 
distributions  with  equal  point  proton  and  neutron  densities  in  both  momen¬ 
tum  and  coordinate  space  formulations  accurately  predict  the  data,  although 
the  two  formulations  require  different  energy  shifts  to  do  so.  This  difference 
is  thought  to  be  a  result  of  the  more  acciuate  nonlocal  representation  of  the 

nuclear  potential  in  the  momentum  space  code.  Additional  spectra  were  ob- 

/ 

tained  for  scattering  at  constant  momentum  transfers  of  .94  and  1.5yfm"^  in 

/ 


order  to  generate  constant  momentum  transfer  excitation  fxmctions.  Use  of 
these  excitation  fimctions,^~|^  ratios,  and  shell  model  DWIA  calculations 
allowed  identification  of  several  excited  states  having  shell-model-like,  sin¬ 
gle  particle-hole,  pure  spin-fiip  excitations.  Shell  mode^  and  collective  model 
DWIA  calculations,  as  well  as  the  q  =  .94  and  1.5yfm''^^xcitation  fimc- 
tions  and  the'  ||^|  ratios  indicate  that  the  other  states  are  generally  well 
represented  by  4  shell  model  description  with  collective  enhancements. 


Chapter  1 


Introduction 

In  1935  Yukawa[YU35]  predicted,  by  a  field  theoretical  treatment,  the 
existence  of  a  particle  of  approximately  100  MeV/c^  mass  as  the  field  quantum 
of  the  nuclear  force.  This  particle,  the  pion  (jt),  was  first  observed  by  Powell 
et.  al.  at  the  University  of  Bristol  in  England  as  a  component  of  cosmic  radi¬ 
ation  using  photo  emulsions.  The  pion  was  first  artificially  produced  in  1948 
by  Gardner  and  Lattes  at  the  University  of  California  Radiation  Laboratory’s 
Synchrocyclotron  facility.  It  very  quickly  became  evident  that  pions  had  many 
advantages  in  their  use  for  the  study  of  nuclei  and  nuclear  interactions  because 
of  their  interaction  via  the  strong  force,  but  there  were  no  facilities  for  pro¬ 
duction  of  pions  on  a  scale  large  enough  to  do  these  studies.  Initial  studies 
therefore  concentrated  on  pion  production  cross  sections  and  basic  properties. 
This  was  the  case  until  the  early  1970’s,  when  the  “Meson  Factories”  at  the 
Los  Alamos  Meson  Physics  Facility  (LAMPF)  in  the  United  States,  the  Swiss 
Institute  for  Nuclear  Research  (SIN)  in  Switzerland,  and  the  Tri-University 
Meson  Facility  (TRIUMF)  in  Canada  came  into  being. 

These  facilities  produced  relatively  high  current  pion  beams,  and  en¬ 
abled  high  resolution  studies  of  nuclear  structure  using  the  pion.  The  experi¬ 
ment  presented  here  is  one  of  hundreds  performed  to  date  at  LAMPF  in  pursuit 
of  a  knowledge  of  nuclear  structure,  and  it  specifically  targets  the  ^®N  nucleus. 
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A  knowledge  of  the  basic  properties  of  the  pion,  along  with  some  of 
the  characteristics  of  its  interactions  with  nucleons  and  nuclei  is  necessary  to 
understanding  its  ad  /antages  as  a  probe  of  nuclear  structure.  The  pion  is  the 
field  quantum  of  the  strong  or  nuclear  force.  Specifically;  it  is  the  mediator  of 
the  longest  range,  attractive  portion  of  this  force,  with  heavier  mesons  mediat¬ 
ing  the  shorter  range  repulsive  forces  of  the  nuclear  core.  It  is  a  pseudoscalar, 
spin-parity  (J^)  =  0“  meson.  Herein  lies  the  first  of  its  advantages  as  a  nu¬ 
clear  probe.  The  fact  that  it  has  spin  0,  as  opposed  to  i  for  other  probes  in 
popular  use,  reduces  the  number  of  individual  initial  to  final  state  scattering 
amplitudes,  so  that  angular  momentum  considerations  lead  to  a  much  simpler 
two-body  scattering  amplitude. 

The  pion  has  3  charge  states,  7r+,  jt®,  and  7r“,  and  so  has  isospin  T=l. 
Because  of  this,  states  of  higher  isospin  are  obtainable  with  pion  scattering  than 
can  be  reached  with  other  probes.  One  of  the  states  that  can  be  reached  using 
the  pion  is  the  A(1232)  or  A(3,3)  resonance,  where  1232  is  the  energy  of  the 
resonance  in  MeV,  and  (3,3)  represents  the  spin  and  isospin  (2J=3,2T=3)  of  the 
state.  The  isospin  identification  of  this  resonance  is  made  based  on  conclusions 
drawn  from  the  assumption  of  isospin  invariance  for  the  strong  interaction  and 
isospin  algebra  performed  for  7r-nucleon  (tt-N)  scattering  through  this  state. 
With  the  (experimentally  supported)  assumption  of  isospin  invariance  of  the 
strong  interaction  (and  ignoring  the  coulomb  interaction)  comes  the  require¬ 
ment  that  the  amplitudes  for  scattering  from  the  nuclear  potential  depend  only 
on  the  system’s  total  isospin[LE73].  So,  for  states  obtained  in  tt-N  scattering 
the  cross  sections  are  only  dependent  on  scattering  amplitudes,  denoted  by  the 
symbol  /,  associated  with  \f„  =  |,  |.  These  amplitudes  are  /i  and  /|, 
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respectively.  Because  of  conservation  of  isospin  and  charge  there  is  nonzero 
amplitude  only  when  isospin  and  z  component  of  isospin  are  both  conserved: 

<r,r,|/Tin7’;>^o  (i.i) 


only  iiT  =  T'  and  =  T^.  Now,  with  the  conventions: 

ir,r,  >  =  |1,  1  >  forTT'*’ 

|1,  0>  for7r“ 

|1,  — 1  >  for  7r“ 

|i,  for  proton  (p) 

^  >  for  neutron  (n) 

one  can  calculate  the  scattering  cross  section: 

o-fT+p  =  I  <  tt+pI/Itt+p  >  f  (1.2) 


and 


<  tt+pI/Itt+p  >  =  <  i,i|  <  1  >  |i,l  > 

3  3  3  3 

^  2’  2  "^  2’  2  ^ 

From  equation  1.1,  /  must  equal  /a,  so  that 


<7*+p  oc  1/3  p 


(1.3) 


(1.4) 


One  may  also  calculate,  by  applying  the  isospin  raising  operator  to  the  T  =  |, 
Tz  =  \  state,  the  T  =  |  components  of  the  7r“p  cross  sections: 


<^7r-p-.7r-p  =  1  <  TT  p|/|  [tT  p  >  p 


(1.5) 
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and 


<^ir-p-*ir<>p  =  I  <  Jr'pl/i  \ir°p  > 


(1.6) 


raising  the  1|»“  f  >  state: 

3  3 

^2~  2^  a|;r-p  > +^|7r'’n  > 


=  “|l.-l>lj.j>+/9|1.0>l5ri> 


H-l> 

1 

h\j2T{iT)  +  1 


(r+  +  r+)  |r(T),r,(;r-)  >  |r(n),r,(n) 


0  >  i  >  +A11,  -1  >  |i,  5  >1 


.1  1 


a^2r(7r)  +  1 

=  > +^k',p> 

so  for  the  amplitudes: 


<7r-p|/3|7r-p>=  Ji/ji 


5^1 


^  -  ir  I  o  ^  t  f 

<’’  PlAI'’‘>=V3-f|V3  =  T^t 


(1.7) 

(1.8) 


For  the  fi  amplitudes,  simply  looking  up  the  Clebsch- Gordon  Coefficients  in¬ 
stead  of  using  isospin  algebra: 


<7r-p|/il7r  p>  = 
<7r-p|/ik°n>  = 

so  for  the  total  cross  sections  come  the  relations: 


(1.9) 

(1.10) 


.ff+p  Ot  1.^2 1 


(1.11) 
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O’lr-p-^ir-p  OC  |g/|  +  g/^P 

_  ,  |2 

<^ir-p-nr«n  OC  I  g  /f  3  I 


(1.12) 

(1.13) 


If  the  A(1232)  resonance  has  T  =  ~  only  the  /a  amplitudes  should  be  involved, 
and  the  result  should  be: 


“^jr+p-nr+p  •  ^it~p-*ir~p  *  ^jr-p-*ff»n  —  9:1:2  (1.14) 

which  is  the  experimentally  observed  result  (see  figure  1.1),  so  that  T  —  ^  con¬ 
tributions  to  this  resonance  are  very  small[LE73].  Similarly,  neutron  scattering 
ratios  are  seen  to  be: 


nr~n  •  ^ir+n-nr+n  •  ^5r+n-»ir*p  —  9:1:2  (1.15) 

These  strong  ratios  for  neutron  and  proton  scattering  are  the 
basis  of  the  third  reason  that  pion  scattering  is  so  attractive  for  use  in  nu¬ 
clear  structure  studies.  An  optical  potential  for  TT-nucleus  (tt-A)  scattering 
can  be  described  in  terms  of  the  tt-N  scattering  amplitudes  using  an  impulse 
approximation[LE80,  LE80A].  This  fact  means  that  the  same  large  ratios 
should  be  present  in  tt-A  inelastic  scattering  that  is  dominated  by  pure  proton 
or  neutron  excitation.  This  has  been  shown  to  be  true  for  pure  proton  and  pure 
neutron  transitions  in  ^^C[DE79]  and  ^^C[H085].  Pion  scattering  therefore  aids 
in  assessment  of  the  extent  to  which  the  shell  model  is  valid,  separation  of  neu¬ 
tron  and  proton  contributions  to  nuclear  excited  states,  study  of  how  important 
collective  enhancements  are  to  shell  model  states,  and  identification  of  purely 
collective  states. 


(qui)  iwxxl£> 
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Figure  1.1:  tt^-N  cross  sections  as  a  function  of  pion  energy. 

The  solid  line  is  the  7r+  —  p  or  tt”  —  n,  the  chain-dashed  line  is  the  total  7r+  —  n 
or  7r“  —  p,  the  dashed  line  is  the  tt"*"  —  n  7r°  —  p  or  tt”  —  p  — >  7r°  —  n,  and  the 
dotted  line  is  the  tt'*'  —  n  or  tt”  —  p  inelastic  cross  section  for  the  production  of 
several  pion  states  or  heavier  particles[KA64]. 
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Siciliano  and  Walker[SI81]  have  shown  that,  assuming  a  single  step 
reaction  mechanism  and  using  the  impulse  approximation,  pion  inelastic  scat¬ 
tering  can  be  used  to  assess  the  relative  contributions  of  spin-flip  (A^  =  1)  and 
non  spin-flip  (A5  =  0)  components  in  transitions  to  the  nuclear  excited  states. 
The  variation  of  cross  section  with  energy  at  a  constant  momentum  transfer, 
go,  near  Qmax  is  given  by: 

^  =  r(£)(n(£,«)cos=((l)  +  «£,«)sin^(«)J  (1.16) 

where  r(E)  is  the  overall  energy  dependence,  n(E,0)  is  the  non  spin-flip  depen¬ 
dence,  and  ^(E,<?)  is  the  spin-flip  dependence.  At  constant  q,  9  must  decrease 
as  E  increases,  and  this  results  in  an  increase  in  cross  section  with  increasing 
energy  (decreasing  0)  for  non  spin-flip  (AS=0)  transitions,  and  the  reverse  for 
spin-flip  (AS=l)transitions. 

This  dissertation  presents  the  application  of  scattering  and  its  listed 
advantages  to  the  study  of  the  nucleus.  The  shell  model  structure  of  this 
nucleus  is  a  very  stable,  filled  p-shell  core  with  a  single  proton  hole  in  the 
outer  Pi  shell,  as  shown  in  figure  1.2.  Considering  this  structure,  one  would 
expect  one  of  the  lowest  energy  states,  that  formed  when  the  proton  hole  in 
the  pi  shell  moves  into  the  p|  shell,  to  be  a  very  strong  candidate  for  a  pure 
proton  transition.  This  is  because,  in  the  extreme  shell  model,  it  involves 
only  the  proton  hole  and  neutrons  cannot  be  involved.  However,  the  work  of 
Macauley  et.  a/.[MA76],  Suzuki[SU76],  and  Horikawa  et.  a/.[H077]  showed 
the  need  for  effective  charge  enhancements  to  cross  section  calculations  for  this 
state,  which  represent  coupling  of  the  shell  model  single  hole  wave  functions 
to  collective  vibrations  of  the  core,  tt*  scattering  data  on  this  single  hole 
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Figure  1.2:  shell  model. 

state,  which  is  a  J"  =  |”  at  6.32  MeV  excitation  energy  are  presented  and 
analyzed  here. 

Large  ratios  have  been  found  in  other  p-shell  nuclei,  but  not 
where  one  would  initially  expect  to  see  them  while  considering  a  shell  model 
representation  like  figure  1.2.  In  three  good  neutron  single  particle  states 
were  studied  using  tt-A  scattering.  In  the  shell  model,  these  states  are  formed 
by  moving  the  only  valence  nucleon,  a  Ipi  neutron,  to  the  Ids,  2s^,  or  Ida 
shells.  If  this  model  were  completely  valid,  a  =  9  ratio  would  be  expected 
for  these  states;  the  observed  values  were  about  2[DE79].  The  reason  cited 
for  this  is  that  the  ground  state  has  strong  mixing  between  the  lp|  and 
Ipi  orbitals.  ^®N,  eis  mentioned  before,  might  be  expected  to  have  a  large 
ratio  for  the  6.32  MeV  state,  but  other  low  lying  states  in  cannot 
be  viewed  as  simple  single  particle  couplings  to  a  stable  closed  shell  core  as 
can  be  done  in  so  the  extreme  ratios  might  not  be  expected.  It  has 
been  proposed,  though,  that  such  ratios  might  be  seen  based  on  coupling  the 
excited  nucleon  to  a  core  which  has  certain  excitation  due  to  its  isospin  state. 
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Specifically,  looking  at  the  ^®N  ground  state  of  figure  1.2,  if  the  Ipi  proton 
is  excited  to  the  2sld  shell,  the  isospin  of  the  holes  left  in  the  core  (Ti,)  is 
1,  while  excitation  of  a  Ipi  neutron  can  give  either  Th  =  0  or  1.  Using  this 
“Th  model”  Oset  and  Strottman[OS82]  predict  that  the  lowest  lying  positive 
parity  states  are  one-proton-particle  two-proton-hole  states,  and  that  the  next 
higher  set  of  states  will  be  one-neutron-particle  one-neutron-one-proton-hole 
states,  separated  in  energy  from  the  first  set  by  4-5  MeV.  With  the  9:1  ratios 
for  and  fjfrj  >  this  model  can  be  tested.  It  must  be  noted  that  even 
if  this  model  is  expected  to  be  valid,  mixing  between  p  shell  orbitals  in  the 
ground  state  and  collective  enhancements  to  the  shell  model  wave  functions 
may  reduce  the  ratios. 

Large  ratios  have  been  found  in  at  higher  energies  and  at 
large  angular  momentum  transfer  values.  In  fact,  the  largest  ratio  was  found 
for  a  state  at  9.5  MeV  which  has  the  largest  angular  momentum  transfer  value 
possible  in  single  particle  excitation  from  the  p  to  sd  shell;  this  is  called  a 
“stretched”  state.  Specifically,  these  transitions  go  from  the  lp|  to  the  Ids 
shell,  with  spin-flip,  giving  a  total  angular  momentum  transfer  of  four,  and 
an  excited  state  J'  value  of  either  or  |^.  Figure  1.3  shows  schematically 
why  this  might  be  expected  in  There  are  two  representations  of  this  state 
in  the  shell  model.  Figure  1.3  a)  shows  a  configuration  where  a  proton  has 
been  excited  from  the  Ipa  to  the  Ipi  shell,  and  a  neutron  has  gone  from  the 
Ipi  to  the  Ids  shell.  Figure  1.3  b)  gives  the  same  isospin  representation,  but 
proceeds  simply  by  the  promotion  of  a  lp|  neutron  to  the  ld|  shell.  Because 
of  the  nature  of  the  tt-A  interaction,  one  step,  direct  interactions  are  expected 
to  dominate.  Therefore  very  little  of  the  left  configuration  is  expected,  and  the 
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Figure  1.3:  stretched  state  configurations. 


dominance  of  the  right  configuration  results  in  an  approximately  1:9  ratio  for 

• 

In  a  similar  situation  might  be  expected,  and  figures  1.4  a)  through 
e)  show  the  possible  configurations.  Figure  1.4  a)  can  be  described  in  a  weak 
coupling  model  as  a  Ids  proton  coupled  to  the  first  2"^  excited  state  of  ^“*0, 
which  has  two  isospin  components,  Th  =  0  or  1,  leading  to  the  possibility  of 
two  proton  dominated  states  here.  The  configurations  in  figures  1.4  b)  through 
e)  are  based  on  coupling  a  s-d  neutron  to  a  core  excited  to  a  2’*'  or 
level.  Because  of  the  greater  closed  shell  stability  (due  to  the  pairing  inter¬ 
action)  of  the  ^‘‘C  core  as  compared  to  the  ^‘‘N,  the  core  states  would  be 
expected  to  lie  at  a  higher  energy.  All  of  the  core  configurations  except  1.4 
b)  are  reached  by  multi-step  processes,  and  would  not  be  expected  to  contain 
much  strength.  Configuration  1.4  b)  would,  however,  be  expected  to  produce 
two  small  (neutron  dominated)  ratio  states  for  the  Th  =  0  and  Th  =  1 
portions  of  this  configuration. 

A  description  of  the  experimental  apparatus  and  procedure  used  to 
obtain  the  ^®N(7r^,7r=^')^®N*  data  is  given  in  Chapter  II.  Data  reduction  and 
analysis,  along  with  some  general  comments  on  the  results,  are  presented  in 
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Figure  1.4:  stretched  state  configurations. 


12 


Chapter  III.  Theoretical  interpretations  of  the  data  are  discussed  in  Chapter 
IV. 


Chapter  2 


Experimental  Apparatus  and  Procedure 


The  data  presented  here  are  the  results  of  an  experiment  performed 
from  7  to  16  January  1983  at  the  Clinton  P.  Anderson  Meson  Physics  Facility 
(LAMPF)  on  Mesita  de  Los  Alamos,  Los  Alamos  National  Laboratory,  New 
Mexico.  LAMPF  is  a  medium  energy  accelerator  and  basic  nuclear  research 
facility,  the  main  purpose  of  which  is  to  advance  the  study  of  nuclear  structure. 
All  measurements  were  performed  using  the  Energetic  Pion  Channel  and  Spec¬ 
trometer,  which  is  a  good  resolution  (~200  keV)  pion  momentum  dispersion 
channel  and  spectrometer  at  LAMPF.  Following  are  general  descriptions  of  the 
accelerator,  and  the  EPICS  channel,  spectrometer,  detectors,  and  electronics, 
and  their  use. 

2.1  Accelerator 

The  LAMPF  accelerator  produces  high  current  beams  of  medium  en¬ 
ergy  protons  (H'*'),  negative  hydrogen  ions  (H“)  and  polarized  negative  hydro¬ 
gen  ions  (P“).  The  H"*"  beam  is  used  to  produce  the  secondary  beams  of  tt* 
which  were  used  in  this  experiment.  There  are  four  basic  portions  of  the  accel¬ 
erator:  the  ion  sources,  the  750  keV  preaccelerators,  the  drift  tube  linac,  and 
the  side  coupled  linac.  These  will  be  discussed  briefly  here  but  more  details 
may  be  found  in  reference  [LI77]. 
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Figure  2.1:  LAMPF  Facility. 


Three  different  ion  sources  are  used  to  produce  the  three  ion  beams, 
two  of  which  can  be  accelerated  at  any  one  time  (!!■*■  and  H“  or  and 

P“).  The  source,  a  Von  Ardenne  Duoplasmatron,  strikes  an  arc  in  a  100 

micron  H2  atmosphere  between  a  cathode  (pulsed  at  120  Hz  to  -135  V  with 
a  pulsewidth  of  500^s)  and  an  intermediate  electrode.  The  resulting  plasma 
is  separated  into  positive  and  negative  components,  with  the  free  electrons 
accelerated  toward  the  the  intermediate  electrode  and  the  zero  potential  anode 
aperture.  and  H"''  ions  accelerate  out  through  this  same  anode  aperture 

because  of  an  external  extraction  electrode,  biased  to  -25kV,  which  produces 

the  required  field  gradient.  H“  and  P"  sources  are  the  same  as  the  H+,  but 
with  attached  charge  exchange  and  RF  polarizing  cells  to  produce  the  negative 
and  polarized  negative  ions  as  necessary. 
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Each  ion  source  is  housed  inside  a  separate  injector  terminal,  all  of 
which  charge  to  750  kV  (+  for  the  H"*"  beam,  -  for  the  others),  to  give  the 
particles  their  initial  boost  to  750  keV.  The  positive  and  negative  ion  sources 
are  pulsed  180®  out  of  phase,  and  so  the  and  H“  are  timed  to  enter  the  linac 
portion  of  the  accelerator  on  the  appropriate  opposite  phases  of  the  accelerating 
cycle.  The  injector  terminals  are  raised  to  750  kV  by  Cockcroft-Walton  type 
high  voltage  electrostatic  generators.  The  first  linac  portion  of  the  accelerator 
is  a  modified  Alvarez  design  drift  tube  type  with  quadrupole  focusing  inside  the 
drift  tubes.  The  RF  sources  for  this  portion  of  the  accelerator  run  at  201.25 
MHz,  and  the  section  accelerates  the  beam  to  100  MeV.  The  second  linac 
portion,  and  last  portion  of  the  accelerator,  is  a  resonantly  side  coupled  cavity 
linac  which  operates  at  800  MHz  and  accelerates  the  beam  to  800  MeV.  At  the 
end  of  the  accelerator  is  the  beam  switchyard,  which  allows  the  H"*"  beam  to 
continue  undeflected  into  area  A,  the  meson  physics  area,  and  switches  the  H" 
or  P“  beam  to  other  areas.  Figure  2.2  shows  the  switchyard  and  experimental 
areas. 


2.2  EPICS 

The  H'*’  beam,  directed  into  area  A,  is  used  to  produce  secondary 
beams  of  pions.  One  of  these  beams  is  that  which  is  apertured  into  the  EPICS 
channel.  It  is  produced  when  the  800  MeV  protons  collide  with  the  spinning, 
water  cooled,  cylindrical  A-1  carbon  target.  The  EPICS  beam  is  taken  off  at 
an  angle  of  35®  from  the  main  proton  beam,  which  trades  pion  beam  intensity 
(maximum  at  0°)  for  spectrometer  angular  coverage,  while  considering  the  re¬ 
quirements  of  the  other  experimental  areas  in  area  A.  The  whole  EPICS  system 


Figure  2.2;  LAMPF  Experimental  Area  and  Switchyard. 
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Figure  2.3:  Energetic  Pion  Channel  and  Spectrometer. 


consists  of  a  momentum  selecting  and  dispersing  channel,  a  scattering  cham¬ 
ber,  and  a  spectrometer  and  its  associated  detectors,  electronics,  and  analysis 
system. 


2.2.1  EPICS  Channel 

The  EPICS  channel  is  described  in  reasonable  detail  in  reference 
[TH70].  The  function  of  the  channel  is  to  select  pions  of  a  certain  charge, 
and  of  a  small  range  of  momentum  values,  and  disperse  the  selected  pions  in 
the  vertical  direction  according  to  their  momentum  values.  The  channel  and  its 
optics  are  shown  in  figure  2.4.  This  vertical  dispersion  (or  analysis)  of  the  beam 
by  momentum  value  allows  for  much  higher  beam  flux  (200  times  higher)  than 
an  essentially  monochromatic  beam  would  provide,  while  retaining  the  resolu¬ 
tion  of  the  beam,  since  momentum  is  a  well  known  function  of  vertical  position 
at  the  target.  Momentum  selection  and  dispersion  in  the  channel  are  done  by 
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2,4:  EPICS  Channel  and  Optics. 
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Table  2.1;  EPICS  system  specifications*. 


CHANNEL 

SPECTROMETER 

Solid  angle 

3.4  msr 

Solid  angle 

«  lOmsr 

Ap/p 

2% 

Ap/p 

14% 

Beam  size  (horiz.) 

8  cm 

Momentum 

100-750  MeV/c 

Beam  size  (vert.) 

20  cm 

Flight  path 

«  12.5m 

Beam  divergence  (horiz.) 

<  lOriuad 

Dispersion 

4  cm/% 

Beam  divergence  (vert.) 

100  mrad 

Energy 

70-300  MeV 

‘Data  taken  from  reference [LA84]. 

four  dipole  bending  magnets  (BM01-BM04)  that  focus  point  to  point  in  the 
vertical  plane,  and  point  to  parallel  in  the  horizontal  plane.  Focussing  (or  trim) 
magnets  (FM01-FM03)  are  multipole  magnets  which  fine  tune  the  beam  in  the 
channel.  Jaws  FJ01-FJ04  provide  adjustments  for  channel  output  beam  mo¬ 
mentum  range,  intensity,  and  divergence.  Channel  output  beam  characteristics 
with  all  jaws  fully  open  are  given  in  table  2.1. 


2.2.2  EPICS  Scattering  Chamber 

The  EPICS  scattering  chamber  is  vacuum  coupled  to  the  channel, 
and  normally  contains  a  target  ladder  system  which  holds  approximately  beam 
size  solid  targets.  Since  a  gas  target  was  used  in  this  experiment,  the  target 
ladder  was  removed  and  replaced  by  a  cylindrical  cooled  gas  target  vessel.  The 
vessel,  shown  in  figure  2.5,  extended  outside  the  limits  of  the  channel  output 
beam  spot  size  in  all  directions.  It  was  produced  by  the  chemical  deposition 
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Figure  2.5:  EPICS  Cooled  Gas  Target  Vessel  and  Systam. 


of  a  25  micron  thickness  of  nickel  from  a  solution  of  nickel-sulfamate  onto  a 
mandrel.  The  mandrel  was  subsequently  etched  away  to  leave  only  the  vessel 
itself.  Analysis  of  the  vessel  showed  it  to  be  99.5%  nickel,  0.5%  cobalt,  with 
negligible  levels  of  lead  and  copper.  Absolute  gas  pressure  in  the  vessel  was 
monitored  regularly  by  means  of  an  absolute  pressure  manometer  attached  to 
the  vessel’s  gas  input  manifold.  Gas  temperatures  in  the  vessel  were  regularly 
monitored,  too,  by  means  of  five  thermocouple  temperature  monitors  placed  in 
various  locations  from  the  top  to  the  bottom  of  the  target  vessel. 

2.2.3  EPICS  Spectrometer 

The  EPICS  Spectrometer  (figure  2.6)  is  vacuum  coupled  to  the  scat¬ 
tering  chamber.  Because  it  is  not  dispersion  matched  to  the  channel,  the  first 
portion  of  the  spectrometer  is  a  quadrupole  triplet  which  images  the  pion  at  the 
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Figure  2.6:  EPICS  Spectrometer. 
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target  onto  a  set  of  chambers  (the  front  chambers).  This  arrangement  is  used 
to  determine  where  (in  vertical  position)  the  pion  came  from,  and  therefore  de¬ 
termine  its  incident  momentum.  The  triplet  focusses  the  pions,  point  to  point, 
onto  the  wire  chambers  with  a  magnification  of  -1  in  the  vertical  direction,  and 
point  to  parallel  in  the  horizontal  direction,  where  the  horizontal  position  at 
the  front  chamber  is  proportional  to  the  scattering  angle  at  the  target.  After 
the  triplet  are  two  60“  bending  magnets,  BM05  and  BM06,  which  provide  for 
analysis  of  the  scattered  pion.  Another  set  of  wire  chambers  (the  rear  cham¬ 
bers)  identifies  the  pion’s  position  and  trajectory  after  analysis  through  these 
magnets.  Knowing  both  front  and  rear  chamber  position  and  trajectory  and 
the  optics  of  BM05  and  BM06  allows  identification  of  the  scattered  pions  scat¬ 
tering  angle  and  change  in  momentum.  EPICS  spectrometer  performance  and 
acceptance  specifications  are  shown  in  table  2.1. 

2.2.4  EPICS  Detector  System 

Parts  of  the  detector  system  were  mentioned  briefly  during  the  spec¬ 
trometer  description,  as  was  necessary  to  a  description  of  its  function.  In  fact, 
most  of  the  detectors  are  physically  associated  with  the  spectrometer,  and  the 
only  ones  that  are  not  are  those  used  in  beam  monitoring  for  the  purpose 
of  normalizing  experimental  runs  to  normalization  and  acceptance  scan  runs. 
These  beam  monitors  include  an  ion  chamber  called  “Beam  on  Target”  (BOT) 
which  is  located  near  the  A-1  target  and  monitors  the  radiation  from  it.  The 
other  monitors  are  called  ICl  and  1ACM02.  ICl  is  another  ionization  chamber, 
located  in  the  scattering  chamber  along  the  pion  beam  axis  downstream  from 
the  target,  which  directly  monitors  the  pion  beam  (and  included  contaminants 


23 


Figure  2.7:  EPICS  Spectrometer  optics  and  Detector  System. 

of  electrons,  muons,  and  protons)  incident  on  the  target.  1ACM02  is  a  toroidal 
charge  integrating  coil  which  monitors  proton  beam  current  upstream  from  the 
A-1  target.  ICl  was  used  as  the  primary  normalization  detector,  but  for  the 
few  runs  where  there  were  problems  with  ICl,  cross  normalization  to  EOT  and 
1ACM02  was  used. 

The  spectrometer  detector  system  is  shown  schematically  in  figure 
2.7.  As  was  described  in  the  spectrometer  section,  the  front  and  rear  wire 
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Figure  2.8:  Kinematics  calculation  flow  chart. 
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Figure  2.9:  EPICS  wire  chamber  wire  plane  representation. 


chambers  provide  for  calculation  of  reaction  kinematics  via  the  flow  chart  of 
flgure  2.8.  Operation  of  these  wire  chambers  (or  Delay-line  readout  drift  cham¬ 
bers)  is  described  in  detail  in  references  [AT81,  M082].  A  basic  description  is 
given  here.  Figure  2.9  is  a  representation  of  one  wire  plane  that  measures  posi¬ 
tion  in  the  x  direction.  The  anode  wires  of  this  plane  are  connected,  at  regular 
intervals,  to  a  2.5ns/cm  delay  line  which  is  biased  to  approximately  2150V. 
Anode  outputs  are  taken  at  each  end  of  the  delay  line,  converted  to  NIM  level 
fast  logic  pulses,  and  are  used  as  stops  for  CAMAC  Time  to  Digital  Conversion 
(TDC)  units.  The  TDC’s  are  started  by  the  trigger  coincidence  logic  signal 
of  the  spectrometer  system,  S2-S3  (S2  and  S3  are  scintillator  detectors  which 
will  be  discussed  later).  Since  S2  and  S3  are  physically  located  downstream 
from  the  wire  chambers,  up  to  200  meters  of  500  coaxial  cable  was  used  to 
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delay  the  NIM  fast  logic  pulses.  Considering  figure  2.9,  and  assuming  that  the 
cabling  for  both  anode  outputs  are  the  same  length,  if  wire  0  fires,  the  time 
registered  on  the  two  associated  TDC’s  will  be  the  same.  The  TDC  time  for 
the  left  output  (Ti)  minus  the  time  for  the  right  output  (T^)  will  be  0.  So  if 
Ti  —  Tr  =  0  =  Tdiff  for  an  event,  it  was  within  ±4mm  of  wire  0.  In  general, 
then 

X  =  ao  +  ai  (Tdiff)  (2.1) 

where  x  is  the  wire  position,  ao  is  a  constant  offset  reference  position,  and  ai 
is  a  constant  involving  the  delay  per  unit  length  of  the  delay  line. 

These  chambers  are  actually  able  to  determine  position  to  ±. 125mm, 
though,  and  this  is  done  by  using  the  fact  that  the  sum  of  the  two  anode  TDC 
signals  is  equal  to  twice  the  drift  time  of  the  ions  in  the  chamber  from  their 
point  of  creation  (the  pion’s  position)  to  the  anode  wire,  plus  a  constant; 

Tf  +  T,  =  2Tdr  +  c  (2.2) 

The  constant  here  is  easily  determined,  since  it  is  just  the  total  delay  along 
both  of  the  anode  lines,  which  is  obviously  a  constant  independent  of  which 
wire  was  triggered.  This  leaves  a  known  drift  time,  which  can  be  converted  to 
a  drift  distance  through  a  calibration  which  only  assumes  uniform  illumination 
of  the  8mm  wide  drift  cells[M082].  The  only  problem  left  now  is  the  decision 
of  whether  to  add  or  subtract  the  drift  distance  from  the  wire  position  x. 
In  this  experiment,  this  decision  was  taken  care  of  by  placing  another  wire 
plane,  with  anode  wires  positioned  half  way  between  those  of  the  first  plane, 
directly  behind  the  first  plane.  These  measurements  give  x  position  at  the  wire 
chambers,  and  a  similar  setup,  with  anode  wires  at  a  90°  angle  to  those  used 
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for  the  X  measurements  give  y  position.  10  cm  downstream  from  the  first  set 
of  wire  planes  in  the  wire  chambers  is  an  identical  set  of  planes  that  give  x  and 
y  positions  once  again,  so  that  positions  x  and  y  are  known,  and  trajectories 
0{=  dx/dz)  and  <f){—  dy/dz)  can  be  calculated. 

The  detectors  presented  so  far  allow  for  normalizations  between  runs 
and  determination  of  kinematics  for  all  events,  but  an  analysis  of  all  events 
wastes  computer  time  and  tape  storage  space.  The  remainder  of  the  detectors 
allow  us  to  classify  the  events  as  to  whether  or  not  they  are  good  pion  events. 
The  first  of  these  are  scintillators  S2  and  S3.  As  was  mentioned  before,  S2-S3 
is  used  as  a  start  trigger  for  the  TDC’s  on  the  front  and  rear  wire  chambers. 
The  actual  overall  hardware  good  event  trigger  which  causes  a  set  of  event 
measurements  (from  all  system  detectors)  to  be  written  to  tape  is  S2*S3*(F1 
or  F2)*(F3  or  F4)  (where  F1-F4  are  front  wire  chamber  signals)  provided  the 
computer  is  not  already  busy.  Between  S2  and  S3  is  a  slab  of  Incite  which  is 
used  to  “range  out”  protons  before  they  reach  S3,  so  that  proton  events  will 
not  be  written  to  tape.  Nonetheless,  some  protons,  as  well  as  a  large  number 
of  electrons  and  muons,  may  make  it  through  both  S2  and  S3.  These  events 
must  then  be  software  rejected.  S2  and  S3  aid  in  the  software  analysis  of  these 
events.  At  the  lower  spectrometer  momenta,  because  of  their  differing  masses, 
electrons,  pions,  muons,  and  protons  can  be  resolved  based  on  their  velocity 
or  time  of  flight  between  the  scintillators.  So,  software  cuts  are  made  based 
on  the  time  of  flight  between  S2  and  S3.  The  pulse  height  from  a  scintillator 
detector  is,  according  to  the  energy  loss  equation[ME66]: 


dE 

dx 


a  +  b  In 
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=  a  +  b 


(2.3) 


where  a,  b,  and  f  are  constants  which  involve  the  charge  of  the  incident  paricle, 
and  E  and  M  are  the  incident  particle’s  total  energy  and  mass,  respectively. 
Since  energy  loss  in  matter  is  dependent  on  the  charge,  mass,  and  energy  of  the 
particle  traversing  it,  cuts  on  S2  and  S3  pulse  heights  (or  the  geometric  mean 
of  their  pulse  heights)  can  be  used  to  similarly  eliminate  non-pion  events. 


At  high  momentum  values  electron,  muon,  and  pion  times  of  flight 

A 

become  unresolvable.  A  Cerenkov  counter,  located  between  S2  and  S3,  can  be 
used  to  eliminate  electron  events,  but  was  not  used  in  this  experiment  because 
electrons  actually  only  made  up  a  small  fraction  of  the  events. 


The  last  detector  system  is  called  the  muon  rejector  (see  figure  2.7), 
and  is  composed  of  a  variable  thickness  aluminum  stack  between  scintillators 
S3  and  S4,  and  carbon  wedges  interspersed  between  scintillators  S4-S9,  along 
with  S4-S9  themselves.  This  system  works  based  on  the  fact  that,  because 
pions  are  lighter  than  muons,  a  pion’s  total  energy  to  mass  ratio  for  a  given 
spectrometer  momentum  value  is  higher  than  a  muon’s,  and  since  energy  loss 
in  matter  goes  as  in  equation  2.3,  the  pion  loses  its  energy  before  the  muon. 
By  placing  the  proper  amount  of  aluminum  and  carbon  between  S3  and  one 
of  S4-S9,  it  can  be  arranged  so  that  only  muons  arrive  at  the  scintillator  of 
interest.  Rejection  of  anything  that  makes  it  to  that  scintillator  is  therefore 
muon  rejection. 
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2.2.5  Data  Acquisition  and  Analysis  System 

Anode  signal  cables  from  the  wire  chambers  are  capacitively  blocked 
to  protect  the  electronics  system  from  their  high  DC  voltage,  and  terminated 
into  lOOfl  at  the  input  of  an  amplifier-discriminator  which  produces  fast  logic 
pulses  that  are  sent  on  to  the  electronics  system  in  the  EPICS  counting  house. 
Light  signals  from  scintillators  S2-S9  are  optically  coupled  into  two  photomul¬ 
tiplier  tubes  each,  whose  signals  are  also  cabled  into  the  counting  house.  A 
schematic  of  the  electronic/logic  system  to  which  these  signal  cables  are  input 
is  presented  in  figure  2.10.  This  schematic  does  not  include  the  elastic  fast 
clear  system  used  in  this  experiment,  which  uses  the  obvious  correlation  be¬ 
tween  crude  x  position  in  the  front  and  rear  chambers  to  reject  9  out  of  10 
elastic  events.  This  system  was  used  on  runs  with  overly  high  elastic  count 
rates  (basically  the  forward  angle  runs)  in  order  to  get  the  inelastics  without 
having  to  run  for  an  unreasonable  period  of  time.  When  the  logic  system  pro¬ 
vides  a  hardware  trigger,  a  PDP  11-45  computer  reads  the  CAMAC  TDC’s 
and  Analog  to  Digital  Converters  (ADC’s,  used  to  read  scalers  like  those  from 
ICl  or  1ACM02).  In  addition  to  this,  the  computer  reads  various  other  scalers 
(like  EVENT  and  EVENT-BUSY).  The  scalers  are  indicators  of  the  system’s 
performance  and  were  monitored  at  regular  intervals  by  the  experimenters.  In¬ 
terface  between  the  CAMAC  modules  and  the  computer  is  provided  by  a  Los 
Alamos  developed  Microprogrammable  Branch  Driver.  The  computer  writes 
the  event  data  and  scalers  to  tape  as  they  are  received,  and  analyzes  them  only 
when  it  is  not  busy  reading  and  writing  them. 
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2.3  Experimental  Detail 

This  section  describes  some  details  specific  to  this  experiment.  The 
target  material  was  cooled  nitrogen  gas  isotopically  enriched  to  99.7%  in  the 
isotope  ^®N.  The  nitrogen  gas  was  maintained  at  a  pressure  of  about  13  psia 
and  a  temperature  of  approximately  88  K.  The  gas  temperature  obtained  was 
achieved  through  the  use  of  a  liquid  helium  coldfinger  in  conjunction  with  a  re¬ 
sistive  heater,  both  in  thermal  contact  with  the  target  vessel.  The  temperature 
quoted  was  an  average  of  that  read  from  the  five  thermocouple  temperature 
monitors  placed  throughout  the  target  vessel.  Maximum  variation  of  any  ther¬ 
mocouple  reading  from  the  average  temperature  was  1.8%.  Target  pressure 
and  temperature  were  logged  every  two  hours.  Nitrogen  gas  temperature  and 
pressure  were  not  constant  throughout  the  experiment  because  the  target  was 
purged  and  refilled  a  number  of  times.  Ratios  of  Pressure  to  Temperature 
(P/T)  for  three  different  series  of  data  runs  were  0.1388  ±0.0004,  0.1466 
±0.0007,  0.1485  ±0.0003  psia/K.  The  P/T  ratios  were  used  in  an  ideal  gas 
law  calculation,  along  with  the  beam  and  target  geometries,  to  get  a  target 
thicknesses. 

Yields  for  7r*-p  elastic  scattering  were  measured  by  filling  the  target 
vessel  with  CH4  (methane)  gas.  The  temperature  of  the  methane  gas  was 
approximately  138  K,  giving  a  P/T  ratio  of  0.0916  ±0.0010  psia/K.  Methane 
data  were  taken  at  5°  steps  from  0iab=25‘’  to  50®,  and  at  10°  steps  from  50° 
to  90°.  These  measurements  provided  absolute  normalization  factors  to  known 
cross  sections  as  well  as  spectrometer  acceptance  normalization  factors  which 
were  necessary  because  of  the  spectrometer-target  geometry,  which  caused  an 
angle-dependent  variation  of  effective  target  thickness  (see  figure  2.11). 


Figure  2.11:  Illustration  of  change  in  effective  target  thickness  with  angle  foi 
EPICS  spectrometer  with  cooled  gas  cylindrical  target  vessel. 
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data  were  taken  for  tt'*'  and  ir~  at  T*  =  164  MeV  for  scattering  angles  from 
25°  to  90°,  in  increments  of  5°.  Additional  data  were  taken  for  120  MeV  7r+, 
at  0iab  =  89.9°  and  49.8°,  and  for  T,+  =260  MeV  at  0i^b  =  48.8°  and  28.5°. 
These  data  were  taken  to  measure  the  energy  dependence  of  the  excitations  at 
constant  q  (momentum  transfer)  values.  ®®Ni  contamination  from  the  target 
vessel  seen  in  the  spectra  were  subtracted  out,  and  the  other  target  vessel 
contaminants  were  found  to  have  a  negligible  effect  on  the  spectra  obtained. 


Chapter  3 


Data  Analysis  and  Results 

3.1  Event  Analysis 

The  objective  of  event  analysis  is  to  determine  the  scattering  angle 
and  energy  loss  (“missing  mast'' )  of  the  pion  in  its  interaction  with  the  nucleus. 
The  pion’s  missing  mass  is  the  difference  between  its  incoming  and  outgoing 
total  energies,  and  is  (approximately)  equal  to  the  excitation  of  the  residual 
nucleus.  The  analysis  is  done  by  the  online  computer  in  a  “may  process” 
mode,  and  when  event  count  rates  are  too  high,  ine  runs  have  to  be  replayed 
later,  offline  in  the  “must  process”  mode,  to  analyze  all  events.  The  analysis 
described  here  applies  to  both  on  and  off  line  analysis  modes. 

Analysis  is  performed  by  the  standard  EPICS  analyzer  “Q”[AM79] 
on  an  event  by  event  basis.  Q  contains  a  main  analysis  subroutine,  PR0C06, 
in  conjunction  with  a  display  file  (EXP703.set,  figure  3.1)  and  an  experiment 
specific  set  of  tests  (which  are  in  the  file  EXP703.tst,  figures  3.2,  3.4,  3.5,  and 
3.8).  These  are  used  to  direct  PR0C06  in  performance  of  appropriate  data 
tests  and  setup  of  data  display  histograms,  dotplots,  and  histogram  storage 
files.  PR0C06  also  computes  chamber  and  scintillator  quantities  (front  and 
rear  chamber  positions  and  trajectories,  for  example).  EXP703.tst  is  composed 
of  3  test  loops,  which  are  performed  in  order  by  PR0C06.  In  between  test 
loops,  PR0C06  calculates  chamber  and  scintillator  quantities,  and  reaction 
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iCi39iai9]EXP703.8eT 
) CYCLE  37  CXP703  NIS  REPLAY 

I 

FILE«C139.3i9]EXP703.  DflP 
FREE  ALL 

Ht>393. 8900.  ttSOO.  OOi  XTOT 
Tt-9i  I  MOD  CHHIRB.  ANO  CHK 

Ha-a34. 8000.  taooo.  SOi  Tt4TT0T 
T2^1  I  «OOD  CHHIRB.  ANO  CHK 

H3-8a9.YOOO.110QO.aOi  YTOT 
T3-Y1  I  0000  CHHIRB.  ANO  CHKN 

H4-23A.  9000. 11000. 40l  PHITOT 
T4-9i  i  OOOD  CHHIRB.  ANO  C:« 

H9-237. 9900. 10900.  Ml  THETACHK 
T9-80  I  0000  CHHIRB 

HA-338. 9900. 10900.  Ml  PHICHECK 
TA-08  I  0000  CHHMB 

H7-1A3. 9990. 13000.901  FLR<1) 
T7-93I  0000  EVENT 

HO-IAA.  9990. 13000.901  FLR(3> 
T8-93I  OOOD  EVENT 
H9«a4A.  9900. 13000.41  HHA88 
T9-109I  OOOD  EVENT.  -1/10 
HtO>34A.  9900. 13000.  4i  HMA68 
TlO-llOt  PION.  -1/10 
Hti-a4A.  9900. 11000.4)  HHA88 
Tll-114i  0000  EVENT.  1/10 
CLEAR  ALL 
ITHAT8  ALL  nXMl 


Figure  3.1:  LAMPF  Experiment  703  Data  Display  File,  Exp703.set. 


kinematics.  A  brief  description  of  this  analysis  sequence  follows.  Loop  1  is 
performed  first,  and  consists  of  tests  1  through  16  in  EXP703.tst,  The  loop  is 
basically  a  set  of  two  windows  or  gates,  called  a  box,  placed  around  the  mean 
energy  loss  in  scintillators  S2  and  S3,  and  around  the  time  of  flight  between  S2 
and  S3.  It  is  called  the  particle  identification,  or  PID,  box.  A  dot  plot  of  S23 
Time  of  Flight  versus  Energy  Loss  (Geometric  Mean  of  S2  and  S3)  with  this 
PID  box  placed  on  it  taken  during  replay  of  the  data  is  shown  in  figure  3.3.  If 
an  event  falls  within  this  box,  the  analyzer  continues  on  the  event  by  comput¬ 
ing  particle  positions  and  trajectories  (xf ,  yf ,  0f,  Xr,  yr,  0r»  andi^r)-  After  these 
calculations,  the  analyzer  goes  to  the  second  loop,  tests  17  through  64.  These 
tests  make  sure  that  all  16  wire  planes  fired,  and  that  all  wire  plane  events 
fall  along  a  straight  line  trajectory.  With  these  tests  passed,  the  analyzer 


Figure  3.3:  Experiment  703  Data  Replay  Particle  Identification  (PID)  Box. 


6AT,l41,SBBB,tB14B, 

»17 

GAT.14E.9SSB.1014B. 

|10 

eAT.14l.999t.iai4B. 

Jt9 

EAT.144,aBaB.lB14B. 

t2S 

SAT.  t4S.9S9i, 10141. 

|21 

6AT.t4E,f990.1B14B, 

>22 

SAT. 14B.B9aB. 10140. 

>23 

BAT. 100.0990. 1014S. 

>24 

SAT.2BI.9i9B, 19140. 

>29 

BAT.2S2.9990. 10140. 

>26 

6AT.2S1.999B.1B14B, 

>27 

6AT.2B4.9990.101«. 

>20 

GAT.205.9990. 10140. 

>29 

SAT,2BB.H9B.1B140. 

>30 

SAT.209.9990. 10140, 

>11 

6AT.210,999B.1B14B, 

>32 

SAT.14F.9926.99S2. 

>33 

SAT. 140.9949.10002. 

>14 

GAT,2r2.9990. 10020. 

>39 

SAT.27S.9990.1BS2S. 

>16 

6AT.274. 9090. 10020. 

>37 

SAT.2r9.999B. 10020. 

>16 

AND. 11.29.16.17,10. 

>19 

AND. 11.29.20.27.20. 

>40 

AND. 11. 17. IS, 19.20. 

|4I 

AND. 11.29.10.11.12. 

>42 

AND. 11.21.22.23.24. 

>41 

AND,4B,41. 

>44 

AND.42.43. 

>49 

AHD.44.45. 

>46 

lgR.U,14, 

>47 

AND.46.4r, 

>40 

AND. IB. 19.20.40,49, 

>49 

AND, 17, 19.28,40.49, 

>80 

AND. 17. 18.20.40,40. 

>8« 

AND  17,16,19.40,49. 

>32 

AND.32.33.24.42.44. 

>53 

AND.21.23.24.42.44, 

>94 

AND.31,3a.24,42,44, 

>99 

ftND.21.22.23.42.44. 

>S6 

ANO.26.27,28.41.49, 

>97 

AN0,2S.27,28.41.4S. 

>58 

AND.29.26.26,41.45, 

>59 

AND.29,36,27,41,4S, 

>sa 

AND.30.31. 32.43,44. 

>61 

AHA. 29,11.12.41.44. 

>62 

RNO,39,30.32.43.44. 

>63 

ftND.29.30.31.43.44, 

>64 

END 

ft  MPim  iCHWER  LOOP  2 
MFim 
MinrH 
NFIM 
sftrmft 
t«nm 
BRFim 
■ft  mnm 
c  Mime 

Mime 
3c  Mime 
«  Manic 
X  Mime 
fic  Mime 
9c  Mme 
■c  Mime 
CNECXSUn  (MMS) 

CtCCKSUn  (QIMD) 

FRONT  MIF  BtF  FLRNS  R 
FRONT  DRF  RtF  FLRHE  C 
REAR  MIF  tIF  PLANE  A 
rear  DRF  DIF  PLANE  C 
DRF  DIFF  T9T»PID  ->  DRF 
FRONT  C  CNRNBERS  OKAPID 
FRONT  ft  CHftHBOC  OKtflD 
REAR  C  CHAHBER8  OK^PID 
REAR  A  a«H«ERS  0K«PID 
FRONT  OK 
REAR  OK 

ALL  EXCEPT  QUAD  OKIf  ID  CHWRB 
QUAD  OK 
■■L  OK 

ALLBUT  lA  CHAieOi  EFFICIENCY  TEST 

ALLRUT  2A 

ALLDUT  3A 

ALLBUT  m 

ALLBUT  SA 

ALLBUT  6A 
ALLBUT  9A 
AUBUT  IBA 
ALLBUT  1C 
AUBUT  2C 
ALLBUT  1C 
ALLBUT  C 
AUBUT  9C 
ALLBUT  6C 
ALLBUT  9C 
ALLBUT  IBC 


Figure  3.4:  Experiment  703  Test  File,  Exp703.tst,  loop  2. 
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calculates  the  kinematical  quantities  incident  momentum, 

scattered  momentum,  and  scattering  angle,  and  from  incident  and  scattered 
momentum  and  scattering  angle,  pion  missing  mass.  Event  position  and  tra¬ 
jectory  angles  at  the  target  are  found  by  tracing  backwards  using  the  position 
and  trajectory  angles  of  the  event  at  the  front  chambers  through  the  known 
quadrupole  optics.  For  example,  since  the  magnification  of  the  quadrupoles  in 
the  X  direction  is  -1, 

xtgt  =  -Xf  (3.1) 

There  are  aberrations  in  the  quadrupole  fields,  and  the  actual  calculations  made 
on  these  quantities  involve  fourth  order  polynomials.  Incident  event  momentum 
(Pj)  is  calculated  according  to  the  following  equation  from  Xtgt  using  the  known 
central  momentum  of  the  channel  (Pc,  determined  by  the  channel  magnetic  field 
settings),  and  the  dispersion  of  the  channel,  Dc  =  (^)^>  see  table  2.1. 

Pi  =  Pc  +  x.„(^)  (3.2) 

Where  Xtgt  is  the  positive  or  negative  event  position  on  the  target  relative 
to  the  position  of  channel  central  momentum  on  the  target.  Scattered  event 
momentum  (Pgc)  is  determined  by  an  analysis  using  the  spectrometer  dipole 
magnets  and  is  therefore  a  function  of  positions  before  and  after  analysis  by 
the  magnets,  and  of  the  spectrometer  central  momentum  and  dispersion.  Since 
the  magnification  of  the  dipoles  in  the  x  direction  is  -1,  the  position  at  the  rear 
chambers  should  be  the  negative  of  that  at  the  front  for  events  which  have  the 
same  scattered  momentum  as  the  spectrometer’s  central  momentum: 


Xr  =  -Xf 


(3.3) 
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or  more  generally  for  all  momenta: 


Pk:  =  ^tpP^  +  Ptp  (3.4) 

to  first  order,  where  D,p  =  and  S^p  =  (xfr  +  Xr)  /D*p.  Finally,  scattering 

angle,  9^  is  given  by 

5.c  =  ^.p  +  V  (3.5) 

where  ^.P  is  the  spectrometer  angle. 

Kinematic  calculations  are  done  by  the  subroutine  CALKIN,  using 
Pi,  Psc)  and  to  calculate  incident  and  scattered  event  total  energies  E;  and 
Ejc,  respectively,  and  residual  nucleus  momentum  using  conservation  of  energy 
and  momentum.  Finally,  CALKIN  uses  conservation  of  energy  in  the  form  of 
the  following  equation: 

E,  +  M„c*  =  E„  +  [(P„c)’  +  (Q  +  M„c"f  ]  ’  (3.6) 

where  Mr  is  the  mass  of  the  target  nucleus,  is  the  momentum  of  the  residual 
nucleus,  to  obtain  a  value  for  Q,  the  nuclear  excitation  energy: 

Q  =  [(QEi  -  E,c  +  M„c2)2  -  (P„c)2]  ^  -  M„c^  (3.7) 

Q  is  also  called  the  missing  mass  because  it  is  approximately  equal  to  the 
energy  lost  by,  or  mass  missing  from,  the  scattered  particle,  for  target  nuclei 
much  more  massive  than  the  particle. 

After  the  analyzer  determines  missing  mass,  it  returns  to  loop  three 
in  the  test  file.  This  loop,  called  background  rejection,  places  cuts  on  the  events 
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Figure  3.5:  Experiment  703  Test  File,  Exp703.tst,  loop  3. 


40 


Figure  3.6:  Experiment  703  Xtgt  Histogram. 

analyzed  to  this  point  based  on  their  Xtgt,  ytgt>  ^tgt)  ^check>  and  (^check  values. 
As  before,  the  “tgt”  subscripted  values  are  position  and  trajectory  values  at 
the  target  after  scattering,  ^check  and  ^check  are  difference  angles  between  0  and 
<f>  measured  at  the  front  chambers  and  9  and  <f)  calculated  at  the  front  chambers 
from  measurements  at  the  rear  chambers.  The  main  purpose  of  cuts  on  <?check 
and  <^check  is  to  get  rid  of  events  where  pions  decay  into  muons  between  the 
front  and  rear  chambers.  Examples  of  histograms  used  to  place  these  cuts  are 
shown  in  figures  3.6  and  3.7.  The  last  loop  in  exp703.tst  produces  the  data 
for  histograms  with  the  various  cuts  or  tests  applied. 


Figure  3.7:  Experiment  703  ^tgt  Histogram. 
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Figure  3.8:  Experiment  703  Test  File,  Exp703.tst,  loop  4, 
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3.2  Data  Reduction 

The  important  results  of  replay  are  missing  mass  histograms  or  spec¬ 
tra  like  the  one  shown  in  figures  3.9  through  3.11.  These  histograms  consist  of 
events  which  have  passed  all  test  loops,  and  which  have  specifically  passed  test 
106  or  110  depending  on  whether  the  feist  elastic  veto  weis  used.  The  histograms 
display  number  of  events  versus  energy,  and  the  events  are  binned  into  energy 
groups  40  keV  wide.  The  spectrometer  has  an  angular  acceptance  of  ±1.5°  and 
all  the  events  from  a  single  angle  setting  were  summed  into  one  histogram,  so 
the  events  are  effectively  binned  into  3°  wide  angle  groups  around  the  nominal 
angular  value.  Number  of  counts  in  (or  area  of)  the  peaks  observed  in  these 
histograms  is  used,  along  with  various  normalizations,  momentum  acceptance 
values,  known  absolute  normalization  cross  sections,  etc.  to  get  differential 
cross  sections  for  excitation  of  the  states. 

Peak  areas  for  the  elastics  were  obtained  using  the  program  LOAF 
[SM78],  while  those  for  the  inelastics  were  extracted  using  the  program  FIT 
[MO].  LOAF  used  a  simple  experimental  reference  peak  shape  to  fit  the  data, 
while  FIT  used  an  experimental  line  shape  folded  with  a  gaussian  of  width 
which  varied  with  excitation  energy,  but  not  with  angle,  for  its  reference  shape. 
Fits  were  performed  separately  on  three  regions  of  the  spectra:  the  elastic  (0 
MeV)  region,  the  region  from  4.5  to  14.2  MeV,  and  that  from  14.2  to  19.5  MeV. 
The  regions  were  fit  in  order  of  increasing  energy,  with  iterative  variation  of 
overall  energy  offset  and  peak  energies,  widths,  and  areas  to  obtain  the  best 
overall  fit.  A  linear  background  was  used  for  all  fits.  Example  fits  for  the  three 
regions  are  shown  in  figures  3.9  through  3.11. 
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Figure  3.9:  The  60®  164  MeV  7r+  Spectrum,  along  with  a  fit  to  its  0  MeV 
region. 
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Figure  3.10:  The  60°  164  MeV  7r+  Spectrum,  along  with  a  fit  to  its  4.5-14.2 
MeV  region. 
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Figure  3.11:  The  60°  164  MeV  tt'*'  Spectrum,  along  with  a  fit  to  its  14.2-19.5 
MeV  region. 
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Peak  areas  along  with  various  correction  factors  are  used  to  get  a 
yield  for  scattering  at  the  various  angles: 


YIELD  =  AREA/(CLT  •  DREFF  •  CHEFF  •  LE  •  SF) 


(3.8) 


where  CLT  is  the  percent  computer  live  time,  or  percent  of  the  time  the  com¬ 
puter  is  not  busy,  and  able  to  read  in  events.  DREFF  is  a  correction  for  drift 
efficiency  and  corrects  for  those  events  which  didn’t  pass  drift  difference  tests 
but  did  fire  all  chambers.  CHEFF  corrects  for  imperfect  chamber  efficiency, 
and  is  a  product  of  the  efficiency  of  the  16  wire  planes: 


16 


CHEFF  =  n 

i=l  L 


All  planes  OK 


(3.9) 


All  but  plane  i  OKj 

LE  is  a  loop  per  event  correction  and  is  basically  a  measure  of  what  fraction 
of  the  written  events  were  analyzed.  Since  all  of  this  data  was  replayed  in  a 
“must  process”  mode,  the  LE  was  1.0.  Finally,  SF  is  a  survival  fraction  factor 
which  corrects  for  the  fact  that  a  certain  percent  (based  on  pion  energy)  of  the 
pions  decay  inside  the  spectrometer  after  scattering: 


SF  =  exp^  p  ’’  ) 


(3.10) 


where  L  is  the  path  length  through  the  spectrometer,  M  is  the  pion’s  rest  mass, 
p  is  the  relativistic  momentum,  and  t  is  its  mean  lifetime. 

To  find  the  differential  cross  section,  for  excitation  to  the  states 
seen  in  the  spectra,  define  the  cross  section,  a,  as  the  probability  that  the 
excitation  will  occur  given  a  certain  flux  of  pions  and  a  certain  target  thickness 
.such  that: 


N  =  <rln 


(3.11) 
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Figure  3.12:  Epics  spectrometer  ^-scan. 


where  N  is  the  number  of  excitation  events,  I  is  the  incident  pion  beam  flux, 
and  n  is  the  number  of  target  nuclei  per  unit  area.  Simple  division  of  cross 
section  by  the  solid  angle  of  the  spectrometer  and  multiplication  by  the  lab 
to  center  of  mass  Jacobian,  G,  gives  differential  cross  section  in  the  center  of 
mass: 


_da__N 
dJlcm  In 


(3.12) 


In  this  experiment,  the  factors  AQ,  I,  and  n  are  not  known  absolutely.  The 
design  solid  angle  of  the  spectrometer  is  10  msr,  but  the  exact  value  is  a  function 
of  where  on  the  spectrometer  focal  plane  (the  rear  chambers)the  event  hits,  or 
equivalently,  it  is  a  function  of  Sgp  (see  equation  3.4).  Because  of  this  the  yield 
as  a  function  of  Sgp  was  measured  for  one  of  the  states  near  a  flat  part  of  its 
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angular  distribution  by  varying  the  spectrometer  fields  to  place  the  state  at 
various  positions  on  the  focal  plane  and  therefore  at  various  5gp  values.  A  plot 
of  this  “^-scan”  is  shown  in  figure  3.12. 

The  ion  chambers  (ICl  and  BOT)  and  proton  beam  integrating  cur¬ 
rent  loop  (1ACM02)  give  quantities  which  are  proportional  to  I,  the  pion 
beam  current  incident  on  the  target.  The  proportionality  is  energy  depen¬ 
dent,  though,  because  the  percent  of  the  beam  accepted  by  the  channel  that 
is  actually  pions  is  energy  dependent.  ICl  is  multiplied  by  a  correction  factor, 
FPI,  which  is  an  energy  dependent  beam  pion  fraction  factor,  to  get  an  en¬ 
ergy  independent  value  proportional  to  I.  This  allows  for  scattering  at  various 
energies,  so  that  excitation  functions  at  constant  momentum  transfer  can  be 
determined. 

Delta  scans  and  pion  correction  fractions  allow  removal  of  focal  plane 
position  and  energy  dependences,  but  they  still  do  not  provide  for  absolute 
knowledge  of  AD  and  I.  Because  of  this,  and  the  fact  that  n  is  not  known 
absolutely  (it  is  a  function  of  scattering  angle  for  the  cooled  gas  target  -  see 
figure  2.11)  the  yield  of  a  reaction  for  which  the  differential  cross  section  is 
well  known  must  be  measured  and  used  for  an  absolute  normalization.  164 
MeV  TT^  scattering  from  methane  gas  at  angles  from  25°  to  50°  in  5°  steps, 
and  from  60°  to  90°  in  10°  steps  provided  this.  Comparison  of  these  measured 
TT^-p  yields  to  those  calculated  using  the  pion-nucleon  phase  shifts  of  Rowe, 
Saloman,  and  Landau[R078]  provides  the  absolute  normalization  factor  Nabs 
used  in  calculation  of  absolute  differential  cross  section: 

d<7  Yield  •  G  •  Nabs 


dDcm  Irel  '  nrel  *  Aflrgi 


(3.13) 
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where  values  subscripted  “rel”  are  experimental  data  values  relative  to  the  CH4 
normalization  data  values. 

The  data  presented  here  include  only  statistical  and  fitting  errors  cal¬ 
culated  by  the  fitting  programs.  Additional  errors  that  should  be  included 
are  in  correction  factors,  peak  fitting,  and  absolute  normalization.  Correc¬ 
tion  factor  errors  include  ±3%  each  for  survival  fraction  uncertainty  due  to 
uncertainty  in  path  length  through  the  spectrometer,  uncertainty  in  chamber 
efficiency  measurements,  and  uncertainty  in  beam  monitoring.  Statistical  er¬ 
rors  in  measurement  of  the  acceptance  scan  were  no  greater  than  2%.  Errors  in 
fitting  the  CH4  norm^llization  data  were  no  greater  than  2.5%,  but  there  was  a 
difference  in  peak  fitting  method  which  introduced  an  additional  error  of  2.5%. 
The  hydrogen  elastic  yields  were  obtained  by  simple  determination  of  the  area 
of  the  peak  minus  the  background.  Actual  experimental  yields  were  obtained 
by  fitting  an  experimental  lineshape  on  top  of  the  background  to  the  peak. 
This  leads  to  a  total  estimated  normalization  error  of  3.5%.  Adding  correction 
factor,  statistical  acceptance  scan,  and  normalization  errors  in  quadrature  gives 
a  total  error,  in  addition  to  the  statistical  experimental  peak  fitting  errors,  of 
7%.  Uncertainty  in  the  relative  normalization  between  7r+  and  7r“,  which  is 
composed  of  all  but  the  acceptance  scan,  normalization  run  peak  fitting,  and 
survival  fraction  errors  (which  are  the  same  for  tt'*'  and  tt"),  is  5.5%. 

3.3  Data 

This  section  includes  a  presentation  of  the  data  as  well  as  some  general 
comments  on  it.  A  more  detailed  comparison  of  the  theory  with  the  data  follows 
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in  chapter  IV.  The  ^®N(7r^,7r*y®N*  data,  in  the  form  of  angle  dependent  cross 
sections,  is  presented  graphically  in  appendix  A,  and  is  tabulated  in  appendix 
B.  The  elastic  distributions  show  the  strongly  forward  peaked,  highly  diffractive 
pattern  characteristic  of  scattering  from  a  black  disk  or  optical  potential.  This 
suggests  the  suitability  of  the  optical  potential  to  a  description  of  the  nuclear 
potential.  Figure  3.13  shows  normalized  spectra  taken  at  a  spectrometer  angle 
of  70“  using  both  It  also  shows  the  difference  spectrum  between  the  7r+  and 
7r~  runs.  Large  peaks  in  the  difference  spectrum  indicate  asymmetries  in  the 
neutron/proton  nature  of  the  state  at  the  energy  of  the  difference  peak.  These 
asymmetries  appear  at  about  5.3,  6.3,  7.6,  10.7,  12.5,  and  17.2  MeV.  They 
might  be  expected  in  light  of  extreme  shell  model  descriptions  of  the  low  lying 
states  as  one  particle  two  hole  positive  parity  states,  or,  for  the  negative  parity 
state  at  6.32  MeV,  as  simple  single  hole  state.  However,  it  is  obvious  that 
such  an  extreme  shell  model  description  is  not  adequate  since,  for  example, 
the  6.32  MeV  single  hole  state  should  be  completely  proton  dominated  if  it 
were.  This  is  because  the  extreme  shell  model  represents  this  state  as  a  ps 
proton  moving  to  fill  the  pi  proton  hole  and  leaving  a  pa  hole  behind.  From 
the  TT^  spectra  and  the  cross  section  ratios  (table  3.1),  it  is  obvious  that 
this  does  not  occur,  since  ,  with  the  dominance  of  the  A(3,3)  resonance,  the 
ratio  should  be  9.  This  is  explained  in  a  model  which  includes  collective 
enhancements  to  the  shell  model  transitions.  At  energies  above  about  10  MeV, 
more  collective  nature  should  be  present,  since  the  excitation  energy  is  above 
the  average  binding  energy  per  nucleon  of  approximately  8  MeV.  Single  particle 
(neutron  or  proton)  excited  states  above  this  energy  should  be  less  likely 
because  concentration  of  this  excitation  energy  on  a  single  particle  should  cause 
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Figure  3.13:  70®  7r+  and  tt  spectra,  along  with  the  diiference  spectrum  between 
the  runs. 
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Table  3.1:  Weighted  average  experimental  cross  section  ratios. 


Energy(MeV) 

5.27 

6.32 

7.16 

7.30 

7.57 

8.31 

8.57 

1.15 

1.78 

0.68 

1.01 

0.62 

1.41 

1.87 

Energy(MeV) 

9.15 

9.76 

9.9 

10.7 

11.3 

11.9 

12.5 

1.73 

0.89 

0.43 

7.05 

1.47 

0.74 

2.57 

Energy(MeV) 

12.9 

13.1 

14.1 

14.4 

14.6 

15.0 

16.5 

) 

1.24 

1.67 

0.94 

0.98 

1.03 

0.99 

1.49 

Energy(MeV) 

16.9 

17.2 

17.6 

18.3 

18.7 

18.9 

1.15 

0.79 

0.85 

1.03 

1.25 

1.25 

it  to  become  unbound.  Figure  3.13  and  table  3.1  show  that  there  are  only  about 
4  significant  ratios  above  10  MeV,  out  of  17  states  identified.  One  state 
at  the  lower  end  of  this  excitation  range  has  large  (  9:1)  ^^4)  ratios,  and  it 
is  reasonable  that  it  would  be  found  here.  Large  ratios  also  appear  at 
12.5,  14.1  and  17.2  MeV,  but  as  might  be  expected  considering  their  higher 
excitation  energy,  the  states  seem  to  be  more  collectively  enhanced  and  have 
smaller,  though  significant,  ratios. 

These  higher  lying  one-particle-two-hole  states  \vefe  described  in  the 
introduction  as  stretched  or  spin-flip  states  in  the  extreme  shell  model.  (See 
figures  1 .4  a)  and  b).)  It  weis  also  pointed  out  that  these  states  could  be  iden¬ 
tified  by  the  shape  of  their  constant  momentum  transfer  excitation  functions. 
These  excitation  functions  for  the  states  identified  as  having  spin-flip  transi¬ 
tions,  and  those  for  the  other  states,  are  presented  and  discussed  in  chapter 
four.  The  spin-flip  states  are  also  high  energy  states  (>10  MeV)  which  have 
large  ratios,  as  should  be  the  case. 

Another  artifact  of  the  stretched  nature  of  these  states  is  their  angu- 
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lar  distributions.  They  have  the  greatest  angular  momentum  transfer  they  can 
with  excitation  of  the  pa  particle  to  the  ps  shell,  AJ=4.  The  slowly  varying, 
almost  flat  shape  of  the  angular  distributions  for  the  states  at  10.7,  12.5,  14.1, 
and  17.2  MeV,  as  compared  to  the  other  angular  distributions,  is  characteristic 
of  AJ=4  transitions.  The  other  transitions  have  more  forward  peaked,  diffrac¬ 
tive  angular  distributions,  characteristic  of  lower  angular  momentum  transfer 
values. 


Chapter  4 


Theoretical  Interpretations  of  the  Data 


Pion-nucleus  (tt-A)  scattering  can  be  related  to  7r-nucleon  (tt-N) 
scattering  through  the  impulse  approximation.  In  the  Distorted  Wave  Im¬ 
pulse  Approximation  (DWIA)  the  assumption  is  that  the  amplitude  for  scat¬ 
tering  between  the  pions  and  the  nucleons  in  the  nucleus  is  the  same  eis  that 
for  free  tt-N  scattering,  and  that  the  effects  of  the  nucleon  being  bound  in 
the  nucleus  are  purely  kinematical.  Details  of  the  DWIA  analysis  have  been 
explored  in  depth  elsewhere  (for  example,  see  reference  [KI84]).  Basically, 
though,  the  tt-A  mvdtiple  scattering  amplitude,  Ta,  can  be  related  to  the 
individual  tt-N  optical  potentials  by[FR53]: 


Ta  =  Va  + - - =  yv(  +  T - -  (41) 

A  +  ‘  |^E-HA-K,-hie  ^  ^ 


Where  E  is  the  collision  energy  in  the  tt-A  center  of  mass  system,  Ha  is 
the  nucleEir  hamiltonian  for  the  A  nucleons,  and  the  Vi’s  are  the  bound  tt-N 
optical  potential  operators.  The  impulse  approximation  involves  replacing 


these  potentieil  operators  for  nucleons  boimd  in  the  nucleus  with  the  well 


known  potential  operators  for  pion  scattering  from  free  nucleons. 


As  was  mentioned  before,  the  feature  of  the  tt-N  interaction  which 
dominates  at  the  experimental  pion  energy  of  164  MeV,  is  the  £=1,  A3/2,3/2 
(1232  MeV)  resonance[K069].  In  particular,  the  introduction  showed  that  at 


55 


56 


this  energy  the  the  pion  is  sensitive  to  the  isospin  nature  of  its  target,  and  so 
a  study  of  the  charge  dependence  on  the  cross  sections  gives  clues  to  the 
neutron/proton  nature  of  the  transitions. 


4.1  Theoretical  Models  for  Elastic  Scattering 

Two  different  DWIA  computer  codes  were  used  in  analysis  of  the 
elastic  scattering  data.  One  of  these,  PIPIT[EI76B,FRA],  describes  the  scat¬ 
tering  through  the  solution  of  a  relativistic  Lippman-Schwinger  equation: 

T(k  k'-  kn)  -U(k  k')+  f  k"\  ko)U{k",  k')  ^  , 

^  ^  ^  ^  y  E{ko)  -  E,{k")  -  E^ik")  +  ie  ^  > 

which  is  in  a  momentum  space  formulation.  Because  of  this  momentum  space 

formulation,  PIPIT  allows  for  easy  inclusion  of  nonlocal  operators  (and  so 

nonlocal  potentials,  i.e.  for  the  spin-orbit  interaction,  etc.),  as  compared 

to  coordinate  space  representations.  A  standard  momentum  space  optical 

potential: 

U(k',  k)  =  {/)p(q)t„p(k',  k;  ko)  +  ^„(q)t,rn(k',  k;  ko)}  (4.3) 

developed  according  to  the  formalism  of  Kerman,  McManus,  and  Thaler 
[KE59]  was  used.  The  off-shell  T-matrix  of  equation  4.2  is  obtained  from 
an  on-shell  matrix  generated  from  the  on-shell  phase  shifts  presented  in  table 
4.1.  These  phase  shifts  are  essentieilly  equivalent  to  those  of  Rowe,  Saloman, 
and  Landau[R078].  Damping  of  the  off-shell  T-matrix  weis  achieved  using  a 
model  which  includes  gaussian  functions  for  the  separable  potentials  which 
make  up  the  off  shell  factors  (<7^(k',  k).  These  factors  have  the  form: 

AU'  „•)  - 


(4.4) 
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Table  4.1:  Phase  shifts*  St  used  in  PIPIT  calculations. 


Ecm  (MeV) 

S3I 

P3I 

P33 

Sil 

Pil 

Pi3 

1186 

-10.6 

-3.2 

38.3 

8.9 

-0.9 

-1.0 

1195 

-11.4 

-3.6 

46.6 

9.0 

-0.5 

-1.3 

1202 

-12.1 

-3.9 

53.6 

9.0 

-0.2 

-1.4 

*The  phase  shifts  presented  here  are  those  which  cover  the  range  of  center 
of  mass  energies  that  are  23,  18,  and  28  MeV  below  that  of  otur  experiment 
(1212  MeV),  corresponding  respectively  to  the  energy  shifted  overall  analysis 
performed  here,  and  to  the  two  variations  of  this  energy  shift  shown  in  figures 
4.3  and  4.5. 

where 

gi(/c)  =  (4.5) 

are  the  gaussian  functions,  and  where  the  «’s  are  the  pion  momenta  in  the 
tt-N  cm  frame,  I,  £,  and  j  are  the  tt-N  system  isospin,  orbital  angular  mo¬ 
mentum,  and  total  angular  momentum,  respectively,  and  the  a^’s  are  the 
gaussiein  damping  factors.  The  labels  of  these  two  factors  in  the  fortran  code 
are  DAMPl  for  ao  and  DAMP2  for  ai,  for  the  ^  =  0  (s-wave)  and  £  =  1  (p- 
wave)  components  of  the  transition,  respectively.  Further,  a  three-parameter 
Fermi  (3pF)  nuclear  density  distribution  was  initially  used  in  these  calcula¬ 
tions,  with  the  peirameters  being  the  3pF  nuclear  charge  density  parameters 
presented  in  [DE87].  Later,  more  accurate  point  proton  and  neutron  nu¬ 
clear  matter  densities  unfolded  for  the  finite  size  of  the  nucleons[FRB]  were 
included  because  of  their  more  accurate  description  of  the  nuclear  matter 
density  distribution.  Calculations  presented  here  are  those  made  with  these 
point  densities.  The  differential  eleistic  cross  section  is  calculated  by  squaring 
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the  scattering  amplitude  in  the  normal  maimer: 


(4.6) 


where 

/(«)  =  /.(«)  +  /*(«)  (4.7) 

and  /c(^)  is  the  overall  ir-nucleus  coulomb  amplitude,  and 

LPIMAX-l 

/*(«)  =  (l/2ik„)  £  (2«+l)e*"(S,-l)P,(cos«).  (4.8) 

/=0 

Here  at  is  the  point  charge  coulomb  phase  shift,  LPIMAX  is  the  number  of 
tt-A  partial  waves  used,  and 

St  =  e^‘*'  =  =  rjte^^^  (4.9) 


where  rjt  represents  the  phase  shifts  for  the  imaginary  or  absorptive  portion  of 
the  scattering  (which  at  164  MeV  is  essentially  nonexistent,  so  that  rjt  —  1), 
and  the  Sf  are  the  real  peirts  of  the  phase  shifts.  The  overlining  simply 
indicates  that  the  phase  shifts  eire  from  coulomb  matched  transition  matrices, 
such  that 


6t  =  ^tan“^ 


-2/>T? 


(4.10) 


in  radians,  where  p  is  a  phase  space  factor.  The  T-matrices  used  here  are  just 
those  generated  in  equations  4. 2-4.4,  and  as  indicated,  are  coulomb  matched. 


The  other  pion  elastic  scattering  code  used  here  was  DWPI  [EI74], 
[EI76A].  DWPI  is  a  coordinate  space  code  which  solves  the  Klein-Gordon 
equation: 

-  (pcf  =  {mcy 


(4.11) 
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where  E  is  the  total  energy,  mo  is  the  pion  rest  mass,  and  p  is  the  momentum 
operator,  p  =  ihv>  These  calculations  used  a  standard  Kisslinger  optical 
potential[KI55]: 

V*{r)  =  (kJbfKr)  -  bf WW  •  V)  (4.12) 

\/here  is  the  pion  momentum,  p(t)  is  the  radial  dependence  of  the  nu¬ 
clear  matter  density  distribution,  and  the  complex  constants  bo  and  bi  are 
determined  according  to  the  equations[AU67]: 

bj  =  F(E„  En)  |(A  -  +  ZaS,|  A"'  (4,13) 

bf  =  F(E„  En)  |(A  -  +  Z(2a-3  +  a-.|  (4.14) 

where  F(E,,EN)is  the  kinematic  factor  relating  the  free  tt-N  phase  shifts  to 
those  for  the  bound  nucleons, 

^2T,2J  =  exp(i52T,2j)sin^2T,2J 

and 


(A  —  Z)  — >  Z  and  Z  (A  —  Z)  for  bp  and  bj 

since  the  potentials  differ  for  7r+  and  tt".  The  phase  shifts  <52T,2J 
are  those  of  Rowe,  Saloman,  and  Landau[R078],  and  the  nuclear  matter 
density  was  agmn  a  3pF  distribution  with  parameters  taken  from  [DE87], 
initially.  The  calculations  presented  here,  though,  used  a  3pF  density  dis¬ 
tribution  with  parameters  adjusted  to  reproduce  the  unfolded  point  density 
distribution  used  in  the  PIPIT  calculation.  This  adjustment  was  accom¬ 
plished  by  using  the  portion  of  PIPIT  which  calculates  the  nuclear  density 
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from  the  3pF  distribution,  and  adjusting  the  parameters  until  the  3pF  dis¬ 
tribution  matched  the  point  distribution.  Finally,  both  elastic  calculations 
employed  a  phenomenologically  motivated  energy  shift  which  caused  the  tt- 
N  interaction  to  be  evaluated  at  an  energy  of  20  to  30  MeV  below  that  of 
the  incident  pions.  Cottingame  and  Holtkamp  have  shown[CO80]  that  such 
an  energy  shift  systematically  improves  DWPI’s  description  of  experimental 
angular  distributions  for  target  nuclei  from  to  ^®*Pb.  This  energy  shift 
has  also  been  found  to  improve  PIPIT  calculations.  The  underlying  physical 
feature  which  is  thought  to  produce  this  effect  is  the  Fermi  momentum  of  the 
struck  nucleon. 


4.2  Theoretical  Models  for  Inelastic  Scattering 

In  the  analysis  of  the  inelastic  data,  two  different  DWIA  codes  were 
used,  ARPIN  and  DWPI.  ARPIN[FRA,LE80]  (Argonne  Pion  Inelastic  Code) 
is  a  momentum  space  DWIA  code  that  uses  microscopic  transition  densi¬ 
ties.  The  tt-N  off-shell  transition  matrix  (T-matrix)  is  generated  by  PIPIT. 
This  T-matrix  along  with  initied  and  final  wave  functions  and  their  associated 
transition  densities  are  input  to  ARPIN,  which  then  calculates  the  differen¬ 
tial  cross  sections.  The  shell  model  calculations  which  produced  these  wave 
functions  and  transition  densities  were  performed  by  D.  John  Millener,  and 
were  essentially  the  same  as  those  performed  in  [MI75]. 

DWPI  was  used  for  both  elastic  and  inelastic  calculations,  and 
macroscopic  transition  densities  were  used  for  the  inelastic  transitions.  In 
this  macroscopic  treatment,  the  radius  parameter  in  the  density  distribution 
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vised  is  made  a  function  of  angle: 

'  =  'ofl+E“A»Yv(f)j  (4-15) 

Here,  is  a  linear  combination  of  creation  and  annihilation  operators 
and  bj^^,  which  create/destroy  a  nuclear  excitation  phonon  of  angular  mo¬ 
mentum  A  and  projection  /z. 

aA„  =  ^aA-'  [bv  +  (-rbLJ  (4.16) 

With  this  deformed  radius  pareimeter  substituted  into  the  3pF  density  distri¬ 
bution,  the  program  continues  on  to  calculate  the  nuclear  excitation.  This  nu¬ 
clear  excitation  is  generated  by  the  portion  of  the  optical  potential  calculated 
from  the  deformed  part  of  the  nuclear  density  distribution.  The  calculations 
presented  here  were  actually  done  with  DWFIT,  which  is  DWPI  modified 
by  R.  Eisenstein  to  use  separate  deformation  parameters,  )0p  and  pn  for  the 
individu2il  proton  eind  neutron  contributions  to  the  density  distribution  and 
optical  potential,  and  also  linked  by  S.  J.  Seestrom  to  the  optimization  code 
MINUIT[JA89],  for  optimization  of  the  fit  of  the  calculated  to  the  experi¬ 
mental  cross  section  by  variation  of  the  parameters  I3p  and  /?„• 

Both  macroscopic  and  microscopic  calculations  were  used  here  for  a 
number  of  reasons.  A  comparison  of  microscopic  calculations  to  experimental 
cross  sections  is  important  because  the  microscopic  model  involves  detailed 
nuclear  structure  calculations  and  this  comparison  provides  feedback  on  our 
knowledge  of  this  structvire.  On  the  other  hand,  some  states  are  found  to 
be  enhanced:  they  are  stronger  than  predicted  by  microscopic  calculations. 
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This  is  explained  by  the  macroscopic  (or  collective)  effective  charge  enhance¬ 
ments  to  the  microscopic  transitions  just  covered.  A  completely  macroscopic 
calculation  like  DWPI  may  also  be  used  to  predict  angular  distributions  for 
very  strong,  collectively  enhanced  states.  Finally,  DWPI  calculates  excitation 
transition  strengths,  which  allows  a  comparison  of  neutron  and  proton  transi¬ 
tion  strengths  and  indicates  the  neutron/proton  nature  of  the  transitions  and 
allows  an  almost  direct  comparison  of  proton  strengths  to  7-ray  transition 
strengths  presently  fotmd  in  the  literature. 

4.3  Theoretical  Results 

4.3.1  Elastic  Scattering 

Calculations  of  the  elastic  scattering  were  performed  with  the  two 
different  codes  described  above.  Input  to  the  two  codes  which  produced 
the  calculated  angular  distributions  shown  in  figure  4.1  are  described  below. 
Both  codes  used  essentieilly  equivalent  ground  state  density  distributions. 
As  mentioned  before,  PIPIT  used  point  proton  and  neutron  distributions 
which  were  unfolded  for  the  finite  size  of  the  nucleons,  and  DWPI  used  a 
3pF  fit  to  this  point  distribution.  The  3pF  parameters  of  this  distribution 
were  C=2.334,  T=0.460,  and  W=0.160,  all  in  fm,  where  C  is  the  density 
distribution  radius,  T  is  the  diffusivity  (or  10  to  90%  thickness),  and  W  is 
the  dip  pEirameter.  Comparison  calculations  using  PIPIT  were  made  in  which 
the  only  parameters  that  differed  between  the  calculations  were  the  3pP  fit 
versus  the  point  densities  distribution,  and  no  discernible  difference  was  seen 
between  the  two.  Energy  shifts  used  in  the  two  codes  were  different,  with 
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the  optimum  shift  for  PIPIT  at  -23  MeV,  and  that  for  DWPI  at  -28  MeV. 
The  shift  of  -28  MeV  for  DWPI  agrees  well  with  the  results  of  Cottingame 
and  Holtkamp,  but  considering  the  association  of  this  shift  with  the  Fermi 
momentum  of  the  nucleons,  it  does  not  agree  with  the  calculated  Fermi  energy 
of  the  nucleons  in  a  nucleus  of  23  MeV[PR75].  The  difference  in  energy  shift 
required  by  the  two  codes  might  be  due  to  the  fact  that  PIPIT  has  a  more 
accurate  nonlocal  description  of  the  potential,  and  therefore  of  the  interaction. 
Finally,  optimum  damping  factors  for  the  tt-N  interaction  range  in  PIPIT 
were  found  to  be  DAMPl  =  3  x  10“®  fm“^  £ind  DAMP2  =  2.5  x  10~®  fm“^ 

Tests  of  the  effect  of  •  hanges  to  elastic  parameters  on  the  ineljistic 
cedculations  were  also  made.  In  these  tests  the  energy  shift  was  varied  in 
both  elastic  calculations,  as  well  as  the  nuclear  force  range  factors  DAMPl 
and  DAMP2  in  PIPIT,  around  the  optimum  values  reported  above.  Dis¬ 
torted  waves  from  these  calculations  were  used  in  inelastic  calculations  of  the 
state  at  5.27  MeV.  Angular  distributions  for  both  elastic  and  inelastic  test 
calculations  me  shown  in  figures  4.2  through  4.5.  These  figures  show  that  for 
both  calculations,  changes  to  the  energy  shift  parameter  produce  much  less 
significant  changes  in  tne  inelastic  distributions  than  in  the  elastics.  There¬ 
fore,  simple  optimization  of  these  parameters  for  the  elastics  and  use  of  the 
distorted  waves  produced  in  the  inelastic  calculations  is  appropriate.  The 
same  is  generally  true  for  the  range  factors  in  PIPIT,  but  going  to  the  larger 
momentum  space  reinge  factors  resulted  in  marked  changes  in  the  elastic  and 
inelastic  distributions.  In  fact,  a  DAMPl  value  of  3  x  lO""*  fm“^  produced  a 
divide  by  zero  error  in  running  the  code.  Care  in  the  selection  of  these  range 
factors  is  recommended  since  it  appears  that  the  code  goes  nonlinear  with 


Figure  4.2:  Effect  of  variation  of  energy  shift  elastic  input  parameter  on 
inelastic  calculations  performed  with  DWPI.  Solid  lines,  dashed  lines,  and 
dotted  lines  are  Egh=23,  28,  and  33  MeV,  respectively. 
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Figure  4.3:  Effect  of  variation  of  energy  shift  elastic  input  parameter  to  PIPIT 
on  inelastic  calculations  performed  with  ARPIN.  Solid  lines,  dashed  lines,  and 
dotted  lines  are  Egh=18,  23,  and  28  MeV,  respectively. 
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Figure  4.4:  Effect  of  variation  of  momentum  space  s  wave  damping  factor  elas¬ 
tic  input  parameter,  DAMPl,  to  PIPIT  on  inelastic  cdculations  performed 
with  ARPIN.  Solid  lines,  dashed  lines,  and  dotted  lines  are  DAMP1=  3x  10"®, 
3  X  10"®,  and  1.5  x  10"^,  respectively. 
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Figure  4.5:  Effect  of  variation  of  momentum  space  p  wave  damping  factor 
elastic  input  parameter,  DAMP2,  to  PIPIT  on  inelastic  calculations  per¬ 
formed  with  ARPIN.  Solid  lines,  dashed  lines,  and  dotted  lines  are  DAMP2= 
2.5  X  10“’^,  2.5  X  10“®,  and  2.5  x  10“®,  respectively. 
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large  values. 

4.3.2  Inelastic  Scattering 

■> 

Figures  4.6  through  4.13  show  the  angular  distributions  for  transi¬ 
tions  from  5.27  up  to  18.9  MeV,  Jilong  with  the  associated  theoretical  calcu¬ 
lations.  Following  are  spin  or  spin-parity  assignments  for  the  states  observed 
here,  with  “possible”  assignments  being  those  made  in  the  forthcoming 
inelastics  paper:  5.27-  (f  ■*■),  6.32-  (f  “),  7.16-  (f  ■^),  7.30-  (f  ■^),  7.57-  (1%  8.31- 
(f^),  8.57-  (§■*■),  9.15-  (§"  and  f+),  9.76-  (f),  9.93-  (f“),  11.29(11.24)-  (pos¬ 
sible  I),  11.88-  (f),  12.94(12.92)-  (f~),  13.15-  (possible  f),  14.38(14.24)- 
d"*^),  14.65-  (possible  1"^),  15.03-  (possible  I'*’),  16.46-  (possible  |),  16.91- 
(possible  |~),  17.58-  (|^),  18.27-  (possible  |),  18.70-  (possible  |)  and  18.91- 
(f+|).  (The  energies  given  here  are  those  at  which  the  states  were  ana¬ 
lyzed,  and  when  there  are  two  different  energies,  the  one  in  parentheses  is 
the  energy  of  the  known  state  with  which  the  state  in  the  present  data  is 
identified.)  With  EPICS’  resolution  at  approximately  200  keV,  the  published 
literature[AJ86]  shows  that  states  of  about  12  MeV  energy  and  above  are 
generally  unresolvable  from  neighboring  states,  and  some  of  them  have  mul¬ 
tiple  spin-peirity  assignments.  Many  of  the  high-lying  states  seen  here  were  so 
weak  that  their  statistics  were  not  good  enough  to  warrant  their  presentation. 

The  cross  section  ratios  of  table  3.1  show  that,  at  energies 

below  about  12  MeV,  asymmetries  in  cross  section  are  present  in  almost  ev¬ 
ery  excitation.  This  indicates  that  these  low  energy  transitions  have  distinct 
proton  or  neutron  characteristics  in  their  transition  densities.  This  in  turn 
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Figure  4.8:  DWPI  macroscopic  model  calculation  results  for  states  with  ex¬ 
citation  energies  from  14.6  to  17.6  MeV. 


Figure  4.10:  Results  for  ARPIN  microscopic  model  calculations  with  no  ef¬ 
fective  charge  enhancements  for  states  with  excitation  energies  from  11.2  to 


15.1  MeV. 
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Figure  4.11:  Results  for  ARPIN  microscopic  model  ceilculations  with  effective 
charge  enhancements  for  states  with  excitation  energies  from  5.27  to  8.31  MeV 
(The  8.31  MeV  calculation  has  been  multiplied  by  100  here). 
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Figure  4.13:  Results  for  ARPIN  microscopic  model  calculations  with  effective 
charge  enhancements  for  states  with  excitation  energies  from  16.4  to  19.0 
MeV.  For  the  state  at  18.9  MeV,  the  dashed  line  is  a  J=  |  calculation,  the 
dotted  line  is  a  J=  ^  calculation,  and  the  solid  line  is  the  sum  of  the  two. 
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suggests  that  these  transitions  are  rather  discrete  shell-model-like  transitions. 
One  also  notices  that  above  this  energy,  the  ratios  are  mostly  approximately 
1.  With  no  preferential  excitation  by  tt'*'  or  tt"  here,  it  appears  that  these 
states  are  probably  much  less  discrete,  that  their  excitation  mechanism  is  ba¬ 
sically  macroscopic.  There  should,  of  course  be  high  lying  shell  model  states, 
but  it  appears  that  because  of  the  high  density  of  states  and  the  spectrom¬ 
eter  resolution  of  approximately  200  keV,  an  identification  of  these  states  is 
difficult.  Nonetheless,  because  single  (or  few)  nucleon  interactions  Eire  most 
likely  in  pion  excitations,  it  would  seem  reasonable  to  assume  that  shell  model 
states  would  be  preferentially  excited.  Some  reasonable  correlation  between 
predicted  high  lying  shell  model  states  and  those  experimentally  observed  was 
achieved.  These  high  energy  spin-parity  assignments  are  questionable, though. 

A  study  of  the  excitation  ftmctions  of  these  states  at  constant  mo¬ 
mentum  transfer  (q)  can,  however,  lead  to  identification  of  high  lying  distinct 
single  component  states.  These  AS=1,  “spin-flip”,  or  “stretched”  excitations 
have  been  identified  in  [SE85].  These  states  were  at  10.7, 12.6, 14.1  and  17.2 
MeV  and  contained  stretched  M4  strength.  Figures  4.14  and  4.15  contain 
plots  of  energy  dependent  cross  sections  at  constant  angle  and  momentum 
transfer  (q  value)  for  the  states  identified  as  having  no  enhanced  spin-flip 
excitation,  while  figure  4.16  shows  these  plots  for  the  M4  states. 

Macroscopic  Model-General  Results  The  macroscopic  model,  DWPI, 
was  used  to  generate  anguleir  distributions  for  excited  states  up  to  an  energy  of 
17.6  MeV.  As  strong  states  are  generally  collectively  enhanced,  their  analysis 
with  a  macroscopic  model  calculation  is  interesting  from  a  qualitative  point 
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Figure  4.15:  Constant  momentum  transfer  excitation  functions  for  states  not 
identified  as  having  “stretched”  transitions  with  excitation  energies  from  13.1 
to  18.3  MeV. 
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of  view.  Such  an  analysis  was  performed  on  these  strong  states,  and  in  fact  on 
all  states  observed.  The  plotted  results  on  states  for  which  DWPI  and  ARPIN 
calculations  produced  reasonable  predictions  are  shown  in  figures  4.6  through 
4.8.  The  basic  result  of  this  qualitative  analysis  is  that  reasonable  calculations 
of  the  strong  states  at  5.27,  6.32,  7.16,  7.57,  9.15,  9.76,  9.9,  11.9,  13.1,  14.4, 
14.6, 15.0, 16.5, 16.9,  and  17.6  MeV  were  made.  Angular  momentum  transfers 
used  in  these  calculations  were  AJ  =  3,  2,  3,  3,  2  and  3,  2,  2,  2,  2,  3,  1,  3,  1, 
2,  and  3  respectively.  These  angular  momentum  transfers  are  consistent  with 
published  spin  and  multipolarity  to  the  ground  state  assignments.  With  the 
calculations  reproducing  the  experimental  cross  sections  well  across  the  whole 
first  diffractive  maximmn,  they  tended  to  underestimate  the  first  minimum 
and  the  second  peak  for  all  but  the  9.76  MeV  state. 

Microscopic  Model-General  Results  Microscopic  calculations  made  for 
states  in  are  shown  in  figures  4.9  through  4.13.  Calculations  other  than 
those  shown  were  made  on  other  high-lying  states,  but  their  correlation  with 
the  data  was  questionable,  so  they  were  omitted  from  this  paper.  ARPIN 
was  used  to  perform  cedculations  with  microscopic  transition  densities  for 
excitation  of  states  in  the  nucleus.  These  calculations  were  made  only 
for  the  positive  parity  states,  and  for  energy  levels  up  to  18.9  MeV.  They 
were  restricted  to  a  Ihu  space  and  produced,  for  a  number  of  the  states 
(especially  the  5/2‘*'),  cross  sections  which  significantly  underestimated  the 
data.  It  was  found  that  “effective  charge”  enhancements  to  the  non-spin 
flip  portions  of  the  transitions  were  needed  to  generate  cross  sections  whose 
magnitudes  were  reasonably  close  to  those  of  the  experimental  cross  sections. 
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These  effective  charges  represented  a  collective  enhancement  to  the  discrete 
AJ=3  transitions.  Because  the  space  is  restricted  to  one  unit  of  angular 
momentum  transfer,  some  higher  L  transition  components  are  omitted  from 
the  wave  function,  or,  from  another  point  of  view,  polarization  of  the  core 
nucleons  by  asymmetries  in  the  valence  occurs,  and  is  not  accoxmted  for  in 
the  extreme  shell  model.  Both  viewpoints  lead  to  the  conclusion  that  this 
model  produces  an  underestimation  of  the  strength  and  cross  section  of  the 
transition.  Isoscalar  and  isovector  effective  charges,  Sq  and  ^i,  are  used  to 
represent  this  polarization  of  the  core  in  the  model,  or  equivalently,  to  make 
up  for  the  restricted  model  space.  These  effective  charges  simply  modify 
the  isoscalar  and  isovector  nucleon  densities[CH87]  according  to  the  following 
equations: 

=  (1  +  (4- 17) 

pUO  =  (l  +  ^i)/)i(r)  (4.18) 

These  enhancements  were  applied  to  the  5/2"^  and  7/2‘*'  states  (the  only 
ones  which  have  LSJ=303  components)  and  were  found  to  increase  the  cross 
sections  as  necessary  to  match  the  data.  Manual  variation  of  isoscalar  and 
isovector  effective  charges  produced  best  values  of  6o=0.72  and  (5i=-0.29  to 
enhance  the  AJ=AL=3,  (non  spin-flip),  AT=0  and  AT=1  portions  of  the 
transitions  to  these  states.  Figures  4.9  and  4.10  only  contain  calculations  for 
I"*"  and  1"^  states  with  no  effective  charge  enhancements,  while  figures  4.11 
through  4.13  contain  the  results  for  enhanced  versions  of  these  calculations 
along  with  those  for  states  which  did  not  require  any  enhancement. 

The  microscopic  calculations,  along  with  effective  charge  enhance- 
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Table  4.2:  Neutron(n)  and  proton(p)  One  Body  Density  Matrix  Elements  for 
calculations  made  on  excitation  of  positive  parity  J'  =  -1^  states  in  with 
a  1-^  basis. 


AJ 

(NLJ).- 

(NLJ)/ 

0p| 

0 

0p| 

0p| 

0d| 

0s| 

Opi 

Opl 

ls| 

1 

Opl 

Odf 

Opf 

ls| 

0p| 

Odf 

Opf 

Odf 

Level  6 
Et=17.68 

n  .00 

-.03 

.47 

.00 

.05 

-.27 

-.04 

.42 

.12 

Eb=18.9 

p  .04 

.02 

.06 

-.07 

.04 

.00 

.00 

-.05 

.16 

Et  and  Eb  are  theoretical  and  experimental  energies  in  MeV. 

ments  as  appropriate,  produced  “good  calculations”  of  the  data  for  the  states 
presented  in  Tables  4.2  through  4.4,  which  are  grouped  by  excited  state  J 
value  and  also  give  the  associated  theoretical  energies  and  excitation  levels. 
A  “good  calculation”  reproduced  the  angular  distribution  and  magnitude  of 
the  data,  and  had  a  theoretical  energy  within  1  MeV  of  the  experimental  en¬ 
ergy.  Some  of  the  states  for  which  imsuccessful  identification  attempts  were 
made  with  the  microscopic  model  Had  no  spin-parity  assignments  in  the  lit¬ 
erature,  or  were  unresolvable  from  nearby  states  with  negative  or  unknown 
parity.  Lack  of  identification  resulted  in  the  conclusion  that  the  strength  in 
the  state  was  probably  from  a  negative  parity  state,  or  in  some  cases,  possibly 
from  a  state  whose  energy  was  too  high  to  be  covered  by  the  set  of  OBDME’s 
used.  Neutron  and  Proton  One  Body  Density  Matrix  Elements  (OBDME’s) 
aie  also  given  in  Tables  4.2  through  4.4. 


85 


Table  4.3:  Neutron(n)  and  proton(p)  One  Body  Density  Matrix  Elements  for 
calculations  made  on  excitation  of  positive  parity  =  I"*"  states  in  with 
a  1-^  basis. 


AJ 

(NLJ).- 

(NLJ), 

OpJ 

id 

opj 

Odf 

1 

Opf  0p|  Opf 
ls|  0d|  0d| 

o  o 

2 

Opi 

0d| 

Opf 

ls| 

Opf 

Odf 

Opf 

Odf 

Level  1 

Et=6.50 

n 

.83 

-.14 

.26  .00  -.23 

.01 

.06 

.13 

.00 

-.04 

Ee=7.30 

P 

.00 

-.09 

.25  .03  -.10 

-.12 

.00 

-.03 

.00 

.00 

Level  2 

Et=8.93 

n 

.02 

-.02 

-.02  -.03  .00 

.27 

.71 

.06 

-.06 

.30 

Ee=8.57 

P 

.00 

.25 

-.02  .05  -.15 

.32 

.00 

-.02 

-.06 

.07 

Level  6 

Et= 13.98 

n 

-.26 

.15 

.49  .15  -.04 

-.29 

-.10 

.06 

.14 

.21 

Ee=14.6 

P 

.00 

-.17 

.37  .00  -.21 

-.22 

.00 

-.05 

.06 

-.08 

Level  9 

Et= 15.67 

Ee=16.5 

Ee=17.6 

n 

-.10 

.01 

.11  .04  -.12 

.00 

-.10 

.19 

-.14 

.46 

Ee=18.3 

P 

.00 

.01 

-.20  -.02  .10 

.01 

.00 

.13 

-.04 

.53 

Level  8 

Et= 14.47 

n 

.20 

.09 

-.11  .03  .05 

-.21 

-.10 

-.15 

-.1 

.44 

Ee=18.7 

P 

.00 

-.14 

-.40  -.01  .17 

-.18 

.00 

-.19 

-.10 

.18 

Level  7 

Et=14.29 

n 

.14 

-.21 

-.33  -.10  .30 

-.20 

.13 

.16 

.10 

.04 

Ee=18.9 

P 

.00 

-.35 

.01  -.07  .16 

-.46 

.00 

.12 

.15 

-.25 

Et  and  Ee  are  theoretical  and  experimental  energies  in  MeV. 
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Table  4.4:  Neutron(n)  and  proton(p)  One  Body  Density  Matrix  Elements  for 
calculations  made  on  excitation  of  positive  parity  J'  =  states  in 

with  a  1-^  basis. 


Tir  _  5+ 

“  2 

AJ 

(NLJ)i 

o  o 

2 

Opi  0p| 
0d|  ls| 

Opf 

Odf 

Opf 

Odf 

Opf 

Odf 

3 

Opf  Opf 
0d|  0d| 

(NLJ)/ 

Level  1 
Et=5.16 

n 

-.03 

.05  .01 

.04 

.00 

-.61 

.26  -.19 

Ee=5.27 

P 

.00 

-.55  -.04 

-.03 

-.06 

-.65 

.21  -.16 

Level  2 
Et=7.42 

n 

-.18 

-.62  -.01 

.02 

-.01 

.31 

-.01  -.25 

Eb=7.16 

P 

.00 

-.21  -.02 

.02 

-.32 

-.25 

.00  .08 

Level  6 
Et=13.50 

n 

.28 

-.05  -.02 

.06 

.00 

-.11 

.02  .45 

Ee=14.4 

P 

.00 

-.14  .22 

-.06 

.26 

-.17 

-.03  .30 

Level  7 
Et=14.48 

n 

-.17 

-.18  -.22 

-.21 

.57 

.06 

-.01  -.17 

Ee=15.0 

P 

.00 

-.08  -.03 

-.05 

.57 

-.10 

.09  -.17 

lir  _  7  + 

-  2 

AJ 

(NLJ).- 

3 

Opf  Opf 
0d|  ls| 

4 

Opf 

Odf 

(NLJ), 

Odf 

Level  1 
Et=6.15 

n 

-.87 

.22  -.20 

-.10 

Ee=7.57 

P 

.00 

.21  -.28 

.00 

Level  2 
Et=11.15 

n 

.07 

.03  .21 

.33 

Ee=11.3 

P 

.00 

-.04  -.40 

-.04 

Ex  and  Eg  are  theoretical  and  experimental  energies  in  MeV. 
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4.4  Discussion 

Following  is  a  discussion  of  results  for  each  state  separately  in  light 
of  the  macroscopic  and  shell  models,  and  in  light  of  published  literature  on 
the  states. 

4.4.1  Low  Lying  Positive  Parity  States 

1/2'*'  states 

The  first  1/2'*'  state  is  identified  in  the  literature[AJ86]  as  being  at 
5.30  MeV,  and  was  unresolved  from  the  5/2'*'  state  at  5.27  MeV.  This  seems 
reasonable  in  that  its  published  strength  is  4  orders  of  magnitude  less  than 
that  of  the  5.27  MeV  state[EN79].  It  has  also  been  identified  as  having  a 
predominant  El  strength  in  transition  to  the  ground  state[EN79].  It  has 
been  suggested[OS82]  that  a  limited  idea  of  the  relative  intensities  of  these 
states  could  be  obtained  from  the  (combined)  angular  distribution,  since  the 
state  should  have  a  maximum  where  the  1/2'*'  has  a  minimum.  A  very 
simple  analysis  based  on  this  idea,  where  the  theoretical  5/2'''  calculations 
are  subtracted  from  the  experimental  data  and  the  result  is  extrapolated  to 
16°,  where  the  peak  1/2'*'  cross  section  should  occur,  gives  a  maximum  l/2‘*' 
Tr"*"  cross  section  of  about  .22mb/sr,  and  tt"  of  about  .26mb/sr,  to  give  a 
ratio  (R)  of  0.85.  These  estimates  of  tt'*'  and  tt"  cross  sections  are 
considerably  higher  than  predicted  by  the  ARPIN  shell  model  calculations 
(on  the  order  of  .03mb/sr)  or  Oset  and  Strottman’s[OS82]  (about  .08mb/sr). 
Finally,  the  R  of  0,85  does  appear  to  agree  with  Oset  and  Strottman’s  quan¬ 
titative  calculations,  but  does  not  support  the  qualitative  zissertion  that,  in 
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a  simple  model  where  this  is  mainly  a  (lpi)“*(2sp,  lp-2h,  Ta=1  state  (the 
Tfc  model),  it  should  be  excited  appreciably  more  in  a  reaction  since  a 
proton  excitation  would  necessarily  yield  Ta=1.  Calculations  by  Lie  and 
Engeland[LI70]  do  not  agree  with  this  simple  model,  but  instead  predict  a 
dominant  3p-4h  strength  of  65%,  and  a  small  lp-2h,  Th=l  strength  of  23%. 
Nonetheless,  ARPIN  microscopic  calculations  agree  with  this  assertion,  show¬ 
ing  the  largest  transition  amplitudes  for  this  component,  (with  AJ=1),  the 
proton  amplitudes  dominating,  and  also  producing  a  theoretical  R  of  1.24. 

The  next  1/2+  state  is  identified[AJ86]  with  an  excitation  energy 
of  8.31  MeV.  This  state  was  present  in  the  data,  but  with  poor  statistics, 
especially  for  the  tt"  cross  sections.  Microscopic  DWIA  calculations  repro¬ 
duced  the  shape  of  the  7r+  experimental  cross  sections  fairly  well,  but  were 
almost  three  orders  of  magnitude  too  low.  Macroscopic  calculations  failed  to 
reproduce  the  shape  at  all.  Table  3.1  shows  a  R  of  1.41,  which  is  at  odds  with 
the  (Th  model)  theory[OS82]  that  this  T/i=0  state  should  be  preferentially 
excited  by  tt".  Lie  and  Engeland  calculate  a  69%  lp-2h  T/,  =  1  contribution 
to  this  state,  in  disagreement  with  the  T/,  model,  but  considering  that  their 
third  1/2+  state  has  a  large  Tfc=0,  one  might  assmne  that  the  second  and 
third  1/2+  states  are  a  T/,=1,0  pair.  Once  again,  though,  the  ARPIN  micro¬ 
scopic  calculation  agrees  with  this  theory  and  shows  a  dominant  (lpp”^(2si) 
(AJ=0)  neutron  amplitude  for  this  state.  One  cannot  say,  because  of  the  poor 
TT"  statistics  and  the  lack  of  agreement  with  the  7r+  data  magnitude,  that  the 
calculation  is  reasonable.  No  other  low  lying  1/2+  states  were  observed,  and 
the  microscopic  calculations  indicate  cross  sections  which  would  be  too  low 
for  them  to  show  up. 


89 


3/2“^  states 

The  first  3/2+  state  is  identified  as  having  an  energy  of  7.30  MeV 
[AJ86],  and  as  having  a  predominant  El  transition  to  the  ground  state[£N79]. 
The  calculations  of  Lie  and  Engeland  show  large  (53%  or  greater)  lp-2h  T;i=0 
strengths  and  small  (6%  or  less)  T/,=l  strengths  for  the  three  lowest  3/2+ 
states,  in  disagreement  with  the  T/,  model.  Calculations  for  this  state  do  show 
its  largest  amplitude  in  an  El  transition,  in  agreement  with  its  identification, 
but  this  amplitude  lies  in  a  Ipi  to  2si  neutron  transition,  which  according 
to  the  simple  T/,  model,  would  indicate  that  it  is  the  Th=0  member  of  one 
of  the  model  pairs.  According  to  the  model,  the  other  member  of  the  pair 
would  have  to  have  a  lower  energy  and  a  large  Ip^  to  2si  proton  transition, 
but  there  is  no  lower  energy  3/2+  state.  However,  the  microscopic  model  does 
not  predict  the  shape  or  amplitude  of  the  7r+  data  well  at  all,  though  it  does 
agree  with  the  tt"  data,  which  has  poor  statistics.  So  the  calculations  are 
suspect,  and  a  simple  T/,  model  might  do  better  here.  The  data  itself  does 
not  really  support  this  model  either,  though,  with  a  R  of  1.01. 

The  second  3/2+  state  is  at  8.57  MeV[AJ86].  The  calculations  agree 
with  the  shape  of  the  7r+  data,  though  the  magnitude  of  the  prediction  is  a 
little  low,  and  once  again,  the  tt"  data  has  poor  statistics.  The  calculations 
have  a  dominant  Ipi  to  Idi  neutron  amplitude,  with  Ipi  to  ld|  neutron 
and  proton  eimplitudes  both  of  about  hedf  the  magnitude  of  the  dominant 
transition.  These  amplitudes  are  also  difficult  to  reconcile  with  the  simple 
T;,  model.  The  experimental  R  here  is  1.87,  which  does  not  agree  with  the 
Th  model.  We  observed  no  other  low  lying  3/2+  states,  and  the  calculations 
indicate  that  they  could  have  been  overshadowed  by  background  or  nearby 
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stronger  states  in  the  spectra. 

5/2'*’  states 

The  state  at  5.27  MeV  has  been  identified  as  the  first  J’^  =  5/2'*‘ 
[AJ86],  and  is  reached  by  a  predominantly  E3  transition  from  the  groimd 
state  [EN79].  In  support  of  the  model  here  are  calculations  by  Lie  and 
Engeland[LI70]  which  indicate  that  there  is  a  large  5/2'*’,  lp-2h,  T/,=l  ampli¬ 
tude  (A=.77)  in  this  state.  Calculations  by  Alford  and  Purser  agree  (A=.93) 
[AL69].  The  experimental  R  of  1.15  is  also  in  support.  Table  4.4  indicates  ap¬ 
proximately  equal  neutron/proton  contributions  in  the  ARPIN  calculations, 
which  predict  the  data  well,  disagreeing  with  the  simple  model. 

The  second  5/2’*'  state  is  identified  at  7.16  MeV[AJ86],  with  a  pre¬ 
dominant  E3  transition(EN79]  to  the  ground  state.  Lie  and  Engeland’s  calcu¬ 
lations  indicate  a  dominant  (80%)  lp-2h  Tfc=0  component  in  this  state,  which 
makes  the  5.27  MeV  state  and  this  one  the  only  pair  on  which  the  simple 
model  agrees  with  their  calculations.  The  experimental  R  of  0.68  also  agrees 
with  the  Th  model.  The  shell  model  calculations  include  a  dominant  M2  Ipx 

to  Idj  transition,  which  in  light  of  the  E3  identification  in  the  literature, 
2 

calls  them  into  question.  They  also  contain  a  somewhat  lower  amplitude  E3 
Ipi  to  Idi  component,  and  with  the  effective  charge  values  used  on  all  of 
the  E3  components  in  the  calculations,  they  are  only  somewhat  low  on  their 
predictions  of  the  tt'*’  and  tt"  data.  The  neutron/proton  amplitudes  for  this 
transition  are  also  approximately  equal,  putting  the  ARPIN  calculations  in 
disagreement  with  the  T/,  model. 

The  third  5/2'*‘  state  is  at  9.16  MeV[AJ86]  and  was  unresolved  from 
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a  3/2“  state  at  9.15  MeV  in  the  data.  The  DWPI  macroscopic  calculations, 
with  a  combination  of  both  J=:3/2  and  5/2  cross  sections,  reproduced  the 
data  well.  The  microscopic  calculations,  on  the  other  hand,  appeared  not  to 
produce  enough  E3  strength  to  predict  the  5/2''‘  portion  of  the  cross  section. 
They  also  indicated  that  the  fourth  5/2'*’  state,  (seen  in  the  literature  at  10.53 
MeV),  should  have  been  strong  enough  to  appeeur  in  the  data,  but  it  didn’t. 

7/2“^  states 

The  7.57  MeV  state  is  identified  as  the  first  7/2+  state[AJ86],  and 
has  a  predominant  E3  strength  to  the  ground  state[EN79].  Lie  and  Engeland 
calculate  an  almost  complete  dominance  (95%)  of  a  lp-2h,  T;i=0  component 
for  this  state,  disagreeing  with  the  T^  model.  The  microscopic  calculations 
don’t  support  the  model,  with  a  dominant  Ip^  to  Idi  neutron  amplitude. 
The  experimental  R  of  0.62  also  disagrees.  There  are  no  other  low  lying  7/2+ 
states,  which  again  calls  into  question  the  simple  T/,  model. 

4.4.2  High  Lying  Positive  Parity  States 

These  states  are  treated  in  order  of  increasing  energy,  and  identifi¬ 
cations  of  are  generally  questionable.  The  identifications  are  made  based 
on  the  presence  of  a  state  of  that  or  imknown  J”  in  the  literature[AJ86]  which 
is  within  the  spectrometer  resolution  energy  of  the  state  of  interest,  a  good 
match  of  the  predicted  angular  dependent  cross  section  to  the  experimental 
cross  section  and,  for  states  of  energy  up  to  16.5  MeV,  a  theoretical  energy 
within  1  MeV  of  the  experimental  energy.  For  states  above  16.5  MeV,  all  of 
which  (except  the  18.91  MeV  state,  which  we  deal  with  separately)  appear  to 
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be  J=3/2  states,  we  identify  them  with  lower  energy  theoretical  states  since 
the  Ihuj  model  runs  out  of  3/2'*'  states  at  16.5  MeV.  Theoretical  energies 
and  OBDME’s  are  included  in  the  tables  presented  for  the  states’  respective 
identified  J”. 

The  first  of  the  high  lying  positive  parity  states  was  seen  experimen¬ 
tally  at  11.3  MeV.  The  state  at  11.29  MeV  is  identified  in  the  literature[AJ86] 
as  a  J”  =  l/2~,  but  the  experimental  angular  distribution  doesn’t  have  the 
shape  of  a  J=l/2  transition.  There  is  a  nearby  state  at  11.24  MeV  which  has 
a  J  assignment  of  3/2  or  greater,  and  since  the  second  level  7/2'^  theoretical 
calculation  at  11.15  MeV  matches  reasonably  well,  we  associate  it  with  this 
state. 

Considering  the  14.4  MeV  state,  there  is  a  14.38  MeV  state  identified 
as  having  J’’=7/2+[AJ86].  None  of  the  7/2+  microscopic  calculations  match 
the  data,  though  the  5/2+  level  6  calculation  does,  and  is  presented  in  Fig. 
4.12.  Macroscopic  model  calculations  with  AJ=3  and  £ui  excited  state  J  of 
7/2  do  match  the  data,  however,  and  these  calculations  are  presented  in  Fig. 
4.7.  Since  the  R  for  this  state  is  0.98,  we  assume  that  it  is  collective  and  rely 
on  the  macroscopic  model  in  saying  it  is  probably  a  J=7/2  state,  with  AJ=3 
from  the  ground  state. 

The  state  at  14.65  MeV  has  an  imidentified  J”.  Both  macroscopic 
and  microscopic  calculations  for  this  state  indicate  that  its  J  is  3/2,  though 
the  two  calculations  have  very  different  shapes.  With  a  R  of  1.03  it  appears 
that  this  state  is  also  collective,  £uid  we  propose  its  description  as  a  J=3/2, 
with  AJ=1  from  the  ground  state,  based  on  the  macroscopic  calculations. 
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The  state  at  15.0  MeV  is  very  strong  in  the  data,  and  both  micro¬ 
scopic  and  macroscopic  calculations  indicate  that  it  is  a  J=5/2  state,  with 
AJ=3  from  the  ground  state.  The  microscopic  calculation  which  predicted 
this  state  was  for  the  seventh  5/2+  level  in  the  model,  and  had  a  theoretical 
energy  of  14.48  MeV.  The  R  of  0.99  indicates  that  it  is  a  highly  collective 
state. 

The  7r“  statistics  on  the  state  at  16.5  MeV  are  poor,  and  we  base  its 
identification  on  the  7r+  data.  Both  calculations  indicate  that  it  is  a  J=3/2, 
with  a  AJ=1  transition  to  the  ground  state. 

The  state  at  17.6  MeV  is  also  a  strong  state  in  tho  data.  There  is  a 
state  at  17.58  MeV  in  the  literature[AJ86]  identified  as  a  3/2+  state,  and  both 
calculations  on  the  17.6  MeV  state  agree  with  this.  Its  R  indicates  that  it  is 
a  collective  state,  and  the  agreement  with  the  magnitude  of  the  microscopic 
calculations  appears  to  be  completely  accidental,  since  the  theoretical  state  is 
predicted  at  an  energy  of  15.67  MeV.  The  microscopic  calculation  was  made 
solely  to  get  a  confirmation  of  the  experimental  since  we  ran  out  of  3/2+ 
levels  and  no  level  Eissignment  was  possible. 

The  states  near  18.3  and  18.7  MeV  have  unknown  We  iden¬ 
tify  them  here  as  possible  J’^=3/2+  states,  since  their  shapes  seemed  to  be 
described  fairly  well  by  microscopic  3/2+  levels  9  and  8  calculations,  respec¬ 
tively.  Both  of  these  states  have  a  R  of  1,25,  indicating  that  they  are  fairly 
collective,  but  the  macroscopic  model  failed  to  predict  their  angular  depen¬ 
dent  cross  section  shapes. 

Finally,  the  state  at  18.91  MeV  is  described  in  the  literature  as  a 
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J’^=(l/2'*'+3/2'*')  state.  Experimental  R’s  once  again  indicate  that  this  state 
is  collective,  but  again,  the  macroscopic  calculations  fail  to  predict  the  shapes. 
Fig.  4.13  shows  microscopic  cedculations  for  the  theoretical  1/2'^  level  6  and 
3/2'*’  level  7  added  to  obtain  a  reasonable  match  to  the  experimental  data. 

4.4.3  Stretched  (M4)  Excitation  States 

Data  and  calculations  taken  from  [SE85]  for  states  which  are  thought 
to  have  significant  strength  from  a  stretched  transition  are  shown  in  figures 
4.17  through  4.20.  This  work  is  also  based  on  the  experiment  703  data. 
As  mentioned  before,  these  states  are  identified  as  stretched  states  based 
on  their  large  (or  small)  ratios,  on  the  shape  of  their  their  angular 

distributions,  and  on  the  (negative)  slope  of  their  excitation  functions  at 
constant  angular  momentum  transfer.  It  must  be  noted  that  the  12.5  and  14.1 
MeV  7r“  distributions  do  not  have  this  characteristic  M4  shape.  This  leads 
one  to  assume  that  these  states  are  not  simply  the  result  of  M4  strength,  but 
also  of  strength  from  lower  multipolarities.  This  assumption  is  supported  by 
the  smaller  than  expected  ratios  found  in  the  experiment.  The  multiple 
multipolairity  assumption  is  also  supported,  in  the  case  of  the  14.1  MeV  state, 
by  its  almost  fiat  excitation  function.  These  facts  led  to  calculations  of  these 
states  which  included  the  addition  of  C2  and  C3  strength.  These  multiple 
multipolarity  calculations  lead  to  more  pronounced  ratios,  as  well  as 
to  a  reconciliation  of  the  angular  distribution  shape  and  excitation  function 
slope  problems.  Considering  only  the  M4  portions  of  these  calculations  results 
in  larger  ratios  of  approximately  9  for  the  10.7  and  12.5  MeV  states, 
and  approximately  2  and  0.7  for  the  14.1  and  17.2  MeV  states,  respectively. 
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Figure  4.18:  tt*  angular  distributions,  along  with  a  M4  calculation  (dotted 
line),  a  C2  calculation  (chain-dash  lines),  and  the  sum  of  these  two  calcula¬ 
tions  (solid  line),  for  164  MeV  pion  scattering  from  to  the  excited  state 
at  12.5  MeV. 
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Figure  4.19:  tt*  angular  distributions,  along  with  a  M4  calculation  (dotted 
line),  a  C3  calculation  (chain-dash  lines),  and  the  sum  of  these  two  calcula¬ 
tions  (solid  line),  for  164  MeV  pion  scattering  from  to  the  excited  state 
at  14.1  MeV. 
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These  ratios  indicate  an  almost  pure  proton  excitation  of  the  10.7  and  12.5 
MeV  states  via  M4  excitation,  proton  dominance  of  the  14.1  MeV  state,  and 
neutron  dominance  of  the  state  at  17.2  MeV. 

The  DWIA  microscopic  calculations  performed  with  ARPIN,  using 
shell  model  wave  function  calculations  which  were  once  again  essentially  the 
same  os  those  presented  in  [MI75],  were  performed  for  comparison  to  the  data, 
and  these  are  shown  in  figures  4.17  through  4.20.  Based  on  this  comparison, 
isoscalar  and  isovector  spectroscopic  amplitudes  were  deduced.  These  ampli¬ 
tudes  are  comparable  to  form  factors  obtained  from  electron  scattering,  and 
this  comparison  was  good  for  all  but  the  state  at  14.1  MeV. 

FVom  a  consideration  of  figure  1.4,  and  its  associated  description  in 
the  text,  two  low  energy  proton  dominated  states  are  expected,  which  have 
a  d|  proton  coupled  to  the  ^‘‘C  2‘*'  Th=0,l  states.  Shell  model  calculations 
using  a  Ifiw  basis  indicate  that  this  configuration  should  actually  result  in  a 
single  I"*”  state  at  11.64  MeV.  The  next  state  is  predicted  by  the  model 
to  be  at  15.66  MeV,  which  is  to  high  an  energy  to  correlate  to  the  12.5  MeV 
experimental  state.  It  is  found  that  when  a  (1  -t-  model  is  considered, 
which  allows  for  3p-4h  configurations,  splitting  of  the  lowest  1"^  level  into  2 
states  occurs.  Additional  evidence  for  the  presence  of  the  3p-4h  strength  here 
comes,  at  least  for  the  10.7  MeV  state,  from  the  strength  seen  for  this  state 
in  the  ^^C(a,  p)^®N  reaction.  Calculations  which  assume  a  simple  two-state 
mixing  of  the  lowest  lp-2h  and  3p-4h  states  result  in  amplitudes  which,  if  one 
assumes  they  add  constructively  for  the  10.7  MeV  state,  produce  cross  sec¬ 
tions  for  excitation  of  the  10.7  and  12.5  MeV  states  consistent  with  those  seen 
experimentally  in  the  ^^C(a,p)^®N  reaction.  (Which  may  mean  nothing  for 
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the  12.5  MeV  state,  since  it  has  an  essentially  0  cross  section  in  this  reaction, 
which  might  mean  that  it  has  0  3p-4h  amplitude).  The  fact  that  neither  the 
14.1  nor  the  17.2  MeV  states  are  seen  strongly  in  the  3  particle  transfer  reac¬ 
tions  may  mean  that  they  are  lp-2h  states.  Electron  scattering  form  factors 
for  the  17.2  MeV  state  quoted  in  [SE85]  indicate  that  this  transition  has  a 
very  strong  isoscalar  excitation.  This  description  is  exactly  that  of  the  second 
1^  state  in  the  shell  model  calculations,  and  also  correlates  with  the  simple 
description  of  this  state  in  a  weak  coupling  model  as  a  d|  neutron  coupled 
to  a  ^‘‘N  2'*'  state  as  in  figtne  1.4b).  This,  plus  the  identification  [A  J86]  of  the 
10.7-12.5  MeV  pair  as  the  lowest  I'*’  states,  indicates  that  the  14.1  MeV  state 
is  J"  =  I'*’. 

Finally,  sum  rule  fraction  calculations  showed  a  smaller  fraction  of 
the  theoretically  predicted  isoscalar  strength  was  detected  in  this  experiment 
than  that  of  the  isovector  strength.  This  is  thought  to  be  the  result  of  the 
greater  fragmentation  of  the  isoscalar  states,  resulting  from  the  fact  that  mese 
states  can  have  either  AT=0  or  AT=1,  as  opposed  to  the  isovector,  which 
can  have  only  AT=1. 

4.4.4  Negative  Parity  States 

At  the  lower  end  of  the  excitation  energies  in  there  is  only  one 
negative  parity  excited  state,  wliich  is  at  6.32  MeV.  This  is  a  one  hole  state, 
where  a  proton  hole  moves  from  the  Ipi  to  the  Ipa  shell,  and  is  the  only  neg¬ 
ative  parity  state  that  can  be  modeled  in  a  l^w  space.  In  the  extreme  shell 
model,  movement  of  this  proton  hole  should  produce  the  only  strength  in  the 
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transition,  and  we  should  therefore  get  a  9:1  ratio  for  R.  The  experimental  ra¬ 
tio  of  1.78  indicates  that  this  transition  is  collectively  enhanced.  Calculations 
with  neutron  and  proton  effective  charges  are  shown  in  Fig.  4.21.  These  cal¬ 
culations  include  the  expected  Opi  — >  Opi,  but  also  contain  a  (LSJ)  =  (202) 
component  scaled  by  neutron  and  proton  polarization  charges  (1  -f  and 
(1  -f  ^p),  where  =  0.3  and  -  1.0.  This  ic,  in  agreement  with  the  work  of 
Macauley  et  al.  [MA76],  Suzuki[SU76],  and  Horikawa  et.  fll.[H077].  Higher 
lying  negative  parity  states  start  at  about  9  MeV  with  transitions  from  the  Ip 
shell  to  the  lfg-2p  shell.  Macroscopic  DWIA  calculations  were  performed  for 
all  states  seen  consistently  in  the  spectra,  and  identified  the  spin-parity  (J’^) 
of  the  low  lying  (less  than  12  MeV)  negative  parity  states  by  energy  in  the 
literature[AJ86].  These  states  are  at  6.32,  9.15,  9.76,  9.93,  and  11.88  MeV. 
Their  assigned  J'  values  are  3/2",  3/2“,  5/2“ ,  3/2“ ,  and  3/2",  respectively. 
Assigned  multipolarities  are  E2  for  the  6.32  and  9.76  MeV  states,  with  which 
the  a '  'ulations  agree.  Further,  the  calculations  for  the  other  low  lying  neg¬ 
ative  parity  states  all  used  spin  transfers  of  AJ=2.  Tenuous  identification  of 
higher  lying  states  at  13.1  and  16.9  MeV  as  3/2“  is  made  in  the  forthcoming 
paper.  This  identification  is  based  on  their  angular  dependent  cross  sec¬ 
tions’  similarity  to  those  of  the  low  lying  negative  parity  states,  and  on  the 
inability  of  the  positive  parity  microscopic  calculations  to  predict  their  shape. 
Other  states  predicted  by  the  macroscopic  DWIA  did  not  correlate  well  with 
the  data. 
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Figure  4.21:  Calculations  along  with  data  for  the  negative  parity  “single  hole 
state  at  6.32  MeV. 


Chapter  5 


Summary 


Experimental  data  on  pion  elastic  and  inelastic  scattering  from 
were  presented.  Microscopic  and  macroscopic  DWIA  calculations  were  also 
presented,  and  were  shown  to  predict  the  data  reasonably  well.  The  mi¬ 
croscopic  calculations  describe  the  J=5/2  and  7/2  states  very  well  (though 
calculations  seemed  lacking  in  their  representation  of  J=3/2  states  as  can  be 
seen  in  their  description  of  the  7.30  and  8.57  MeV  states).  Nonetheless,  the 
shell  model  predicts  some  states  that  are  more  neutron  or  proton  dominated 
than  we  observe.  Effective  charge  enhancements  applied  here  make  the  shell 
model  calculations  of  <7,r+  antf  <3^^-  more  equal,  and  allow  the  predictions  to 
closely  match  the  experimental  data.  These  isoscalar  and  isovector  effective 
charges  were  found  to  be  ^o=0.72  and  ^j=-0.29,  respectively,  to  enhance  the 
AJ=AL=3,  (non  spin-flip),  AT=0  and  AT=1  portions  of  the  transitions  to 
these  states.  The  complementary  use  of  ratio  comparison,  shell  model 
descriptions  of  the  nucleus,  collective  enhancements  to  the  shell  model  states, 
and  purely  macroscopic  calculations  allow  us  to  determine  much  about  the 
structure  of  the  nucleus. 

Some  tenuous  new  J'  identifications  were  made,  and  are  included 
in  the  following  list  of  states  seen  in  the  experimental  spectrum:  0.00-  (|~), 
5.27-  (§■*■),  6.32-  d"),  7.16-  (f^),  7.30-  (f^),  7.57-  (f^),  8.31-  (i"^),  8.57- 
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(§+),  9.15-  (f  and  f ),  9.76-  (f),  9.93-  (f),  10.69-(f+),  11.29(11.24)- 
(possible  f),  11.88-  (f),  12.49-(f^),  12.94(12.92)-  (f),  13.15-  (possible  f), 
14.09-(possible  1"^),  14.38(14.24)-  (|''’),  14.65-  (possible  f"^),  15.03-  (possible 
I"*"),  16.46-  (possible  |),  16.91-  (possible  |"),  17.15-(possible  I"*"),  17.58-  (f"^), 
18.27-  (possible  §),  18.70-  (possible  f )  and  18.91-  (§+^).  (The  energies  given 
here  are  those  at  which  the  states  were  analyzed,  and  when  there  are  two 
different  energies,  the  one  in  pjirentheses  is  the  energy  of  the  known  state 
with  which  the  state  in  the  present  data  is  identified.)  An  evaluation  of  the 
simple  model  of  low  lying  states  being  split  into  Tfc=l,0  pairs[OS82]  is  also 
made.  It  appears  that  the  T;,  model  is  too  simple,  that  collective  effects  or 
other  more  involved  shell  model  effects  complicate  the  picture  too  much  for 
the  model  to  be  of  use  at  all,  since  its  only  success  is  with  the  lowest  energy 
excited  state,  a  5/2’*'  state,  and  the  next  highest  energy  5/2+  state. 

The  ratios,  constant  momentum  transfer  excitation  functions, 
and  angular  distribution  shapes  were  instriunental  in  the  identification  of 
transitions  which  contained  significant  M4  strength.  These  states  were  found 
at  10.7,  12.5,  14.1,  and  17.2  MeV.  Spectroscopic  amplitudes  were  extracted 
for  these  states,  and  they  were  in  reasonable  agreement  with  the  amplitudes 
extracted  from  electron  scattering.  Shell  model  calculations  performed  with 
a  Ihu  basis  produced  only  one  low  energy  I'*’  state.  A  (l-{-3)^u;  basis  was 
required  to  reproduce  the  two  experimentally  observed  states  at  low  energy. 
The  need  for  this  (1-f  3)ftw  basis  indicates  that  3p-4h  components  are  present 
in  these  states,  and  it  is  thought  that  the  lowest  lp-2h  and  3p-4h  states  mix 
to  form  the  experimentally  observed  states  at  10.7  and  12.5  MeV. 

Finally,  both  macroscopic  and  microscopic  calculations  for  the  single 
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hole  negative  parity  state  at  6.32  MeV  were  performed.  The  ratio  was 
1.78,  contrary  to  the  extreme  shell  model  predictions  that  this  state  should 
be  completely  proton  dominated  and  have  a  ratio  of  9.  It  was  found 
that  neutron  and  proton  collective  effective  charge  enhancements  of  (1  +  ^n) 
and  (1  +  5p),  where  =  0.3  and  8^  =  1.0  were  necessary  to  appropriately 
enhance  the  angular  distributions.  For  the  rest  of  the  negative  parity  states, 
only  macroscopic  calculations  were  performed,  and  these  mostly  reinforced 
the  identifications  found  presently  in  the  literature,  although  some  tenuous 
new  identifications  were  made. 


Appendix  A 


Experimental  Data  Plots 
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Appendix  B 


Tabulated  Data 
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0.00  MeV  (Elastic)  tabulated  data 


error 

(deg) 

(mb/sr) 

(mb/sr) 

25.52 

0.1789E+03 

0.2843E+01 

30.62 

0.8647E+02 

0.3274E+00 

35.71 

0.3234E+02 

0.1121E+00 

40.79 

0.9625E+01 

0.4958E-01 

45.87 

0.1082E+01 

0.2168E-01 

50.94 

0.3421E+00 

0.8070E-02 

56.01 

0.1414E+01 

0.4652E-01 

61.07 

0.2160E-f01 

0.2063E-01 

66.11 

0.2024E+01 

0.3033E-01 

71.15 

0.1458E+01 

0.1543E-01 

76.19 

0.7613E+00 

0.2260E-01 

81.21 

0.2671E+00 

0.6673E-02 

86.22 

0.8064E-01 

0.5922E-02 

91.22 

0.7055E-01 

0.2468E-02 

1T~ 

25.52 

0.1858E+03 

0.1805E+01 

30.62 

0.8261E+02 

0.2860E+0^t 

35.71 

0.2907E-I-02 

0.1602E+00 

40.79 

0.6846E+01 

0.7764E-01 

45.87 

0.5873E+00 

0.2178E-01 

50.94 

0.9418E+00 

0.1671E-01 

56.01 

0.2252E+01 

0.8038E-01 

61.07 

0.2847E+01 

0.4225E-01 

66.11 

0.2424E+01 

0.8978E-01 

71.15 

0.1464E+01 

0.1666E-01 

76.19 

0.7565E+00 

0.4291E-01 

81.21 

0.6014E-01 

0.6127E-02 

86.22 

0.1145E+00 

0.1824E-01 

91.22 

0.1178E+00 

0.6158E-02 

5.27  MeV  tabulated  data 


^cm 

az. 

6m 

error 

(deg) 

7r+ 

25,53 

(mb/sr) 

(mb/sr) 

0.3464E+00 

0.2091E-01 

30.62 

0.4312E+00 

0.1487E-01 

35.71 

0.5264E+00 

0.1422E-01 

40.80 

0.6136E+00 

0.1182E-01 

45.88 

0.6094E+00 

0.1522E-01 

50.9' 

0.5469E+00 

0.9067E-02 

56.02 

0.3805E+00 

0.2321E-01 

61.08 

0.2404E+00 

0.6566E-02 

66.12 

0.1232E+00 

0.7366E-02 

71.17 

0.6969E-01 

0.3271E-02 

76.20 

0.5621E-01 

0.6416E-02 

81.22 

0.6975E-01 

0.3237E-02 

86.23 

0.9747E.01 

0.6048E-02 

91.24 

0.9868E-01 

0.2817E-02 

TC~ 

25.53 

0.3821E+00 

0.3039E-01 

30.62 

0.4435E+00 

0.2064E-01 

35.71 

0.4861E+00 

0.2122E-01 

40.80 

0.5983E+00 

0.2170E-01 

45.88 

0.5766E+00 

0.2002E-01 

50.95 

0.5008E+00 

0.1122E-01 

56.02 

0.3703E+00 

0.3242E-01 

61.08 

0.2396E+00 

0.1196E-01 

66.12 

0.8082E-01 

0.1852E-01 

71.17 

0.4385E-01 

0.2917E-02 

76.20 

0.3687E-01 

0.1205E-01 

81.22 

0.4342E-01 

0.5574E-02 

86.23 

0.4272E-01 

0.1433E-01 

91.24 

0.6651E-01 

0.4605E-02 

6.32  MeV  tabulated  data 


0cm 

dg 

error 

(deg) 

i— + 

(mb/sr) 

(mb/sr) 

TT 

25.53 

0.6591E+00 

0.3161E-01 

30.62 

0.8950E+00 

0.1982E-01 

35.71 

0.8818E+00 

0.1742E-01 

40.80 

0.7543E+00 

0.1250E-01 

45.88 

0.5182E+00 

0.1379E-01 

50.95 

0.3362E+00 

0.7062E-02 

56.02 

0.1859E+00 

0.1622E-01 

61.08 

0.7375E-01 

0.3671E-02 

66.12 

0.3965E-01 

0.4301E-02 

71.17 

0.4128E-01 

0.2516E-02 

76.20 

0.6181E-01 

0.6440E-02 

81.22 

0.5786E-01 

0.2958E-02 

86.23 

0.5256E-01 

0.4622E-02 

91.24 

0.3974E-01 

0.1805E-02 

7r“ 

25.53 

0.4299E+00 

0.3783E-01 

30.62 

0.5851E+00 

0.2199E-01 

35.71 

0.5703E+00 

0.2220E-01 

40.80 

0.4421E+00 

0.1895E-01 

45.88 

0.3188E+00 

0.1501E-01 

50.95 

0.1799E+00 

0.6900E-02 

56.02 

0.9589E-01 

0.1776E-01 

61.08 

0.3053E-01 

0.4744E-02 

66.12 

0.9123E-02 

0.1041E-01 

71.17 

0.2066E-01 

0.2129E-02 

76.20 

0.2020E-01 

0.9353E-02 

81.22 

0.3651E-01 

0.5044E-02 

86.23 

0.1431E-01 

0.9935E-02 

91.24 

0.1515E-01 

0.2332E-02 

7.16  MeV  TT*  tabulated  data 


TT^ 


^cm 

error 

(deg) 

(mb/sr) 

(mb/sr) 

25.53 

0.6042E-01 

0.1279E-01 

30.62 

0.3881E-01 

0.7342E-02 

35.71 

0.3014E-01 

0.5775E-02 

40.80 

0.3105E-01 

0.4398E-02 

45.88 

0.3132E-01 

0.5567E-02 

50.95 

0.2094E-01 

0.2725E-02 

56.02 

0.1649E-01 

0.7537E-02 

61.08 

0.1192E-01 

0.2352E-02 

66.12 

0.6549E-02 

0.3032E-02 

71.17 

0.4138E-02 

0.1047E-02 

76.20 

0.1775E-02 

0.2613E-02 

81.22 

0.3703E-02 

0.1291E-02 

91.24 

0.4840E-02 

0.9693E-03 

TT 


25.53 

30.62 

35.71 

40.80 

45.88 

50.95 

56.02 

61.08 

71.17 

81.22 

91.24 


0.6370E-01 

0.9974E-01 

0.6188E-01 

0.4977E-01 

0.6202E-01 

0.3455E-01 

0.7120E-02 

0.1828E-01 

0.3229E-02 

0.1615E-01 

0.8862E-02 


0.2465E-01 

0.1198E-01 

0.1056E-01 

0.9150E-02 

0.9275E-02 

0.4884E-02 

0.1043E-01 

0.5130E-02 

0.1137E-02 

0.3844E-02 

0.2271E-02 


7.30  MeV  TT^  tabulated  data 


^cm 

d<r 

du; 

error 

(deg) 

TT*^ 

(mb/sr) 

(mb/sr) 

TT 

25.53 

0.4357E-01 

0.1345E-01 

35.71 

0.1353E-01 

0.5829E-02 

40.80 

0.1592E-01 

0.4838E-02 

45.88 

0.2547E-01 

0.6196E-02 

50.95 

0.1846E-01 

0.3146E-02 

56.02 

0.1056E-01 

0.7879E-02 

61.08 

0.2815E-02 

0.2702E-02 

66.12 

0.6747E-02 

0.3715E-02 

76.20 

0.1089E-01 

0.3362E-02 

81.22 

0.8841E-02 

0.1601E-02 

86.23 

0.1584E-01 

0.3207E-02 

91.24 

0.1168E-01 

0.1324E-02 

7r“ 

25.53 

0.1690E+00 

0.2867E-01 

35.71 

0.1027E-01 

0.1099E-01 

45.88 

0.1658E-01 

0.1052E-01 

50.95 

0.2990E-01 

0.5997E-02 

56.02 

0.3447E-01 

0.1517E-01 

61.08 

0.1088E-01 

0.6344E-02 

76.20 

0.4581E-02 

0.8287E-02 

81.22 

0.1715E-02 

0.3007E-02 

86.23 

0.2553E-01 

0.1196E-01 

91.24 

0.1327E-01 

0.2757E-02 

7.57  MeV  tabulated  data 


0cm 

Hz 

du> 

error 

(<leg) 

Tr"*" 

25.53 

(mb/sr) 

(mb/sr) 

0.2142E+00 

0.1652E-01 

30.62 

0.2715E+00 

0.1182E-01 

35.71 

0.3765E+00 

0.1177E-01 

40.80 

0.4560E+00 

0.1008E-01 

45.88 

0.4433E+00 

0.1282E-01 

50.95 

0.4185E+00 

0.7533E-02 

56.02 

0.3138E+00 

0.2060E-01 

61.08 

0.1974E+00 

0.5787E-02 

66.12 

0.8891E-01 

0.6323E-02 

71.17 

0.4255E-01 

0.2493E-02 

76.20 

0.2487E-01 

0.4436E-02 

81.22 

0.4207E-01 

0.2660E-02 

86.23 

0.4672E-01 

0.4794E-02 

91.24 

0.7462E-01 

0.2449E-02 

7r~ 

25.53 

0.4125E+00 

0.3077E-01 

30.62 

0.6178E+00 

0.2214E-01 

35.71 

0.7393E+00 

0.2394E-01 

40.80 

0.7763E+00 

0.2399E-01 

45.88 

0.8022E+00 

0.2316E-01 

50.95 

0.6752E+00 

0.1261E-01 

56.02 

0.5230E+00 

0.3583E-01 

61.08 

0.2848E+00 

0.1309E-01 

66.12 

0.1789E+00 

0.2514E-01 

71.17 

0.6348E-01 

0.3383E-02 

76.20 

0.5616E-01 

0.1347E-01 

81.22 

0.6998E-01 

0.6538E-02 

86.23 

0.5419E-01 

0.1451E-01 

91.24 

0.8085E-01 

0.5193E-02 

8.31  MeV  TT^  tabulated  data 


^cm 

dff 

du 

error 

(deg) 

TT'*' 

25.53 

(mb/sr) 

(mb/sr) 

0.1930E-01 

0.1044E-01 

30.62 

0.6376E-02 

0.5882E-02 

45.88 

0.1005E-01 

0.3510E-02 

50.95 

0.3939E-02 

0.1737E-02 

56.02 

0.1152E-01 

0.6188E-02 

61.08 

0.4903E-02 

0.1431E-02 

66.12 

0.4216E-02 

0.1916E-02 

71.17 

0.1551E-02 

0.8509E-03 

76.20 

0.1105E-02 

0.1997E-02 

81.22 

0.3155E-02 

0.9353E-03 

86.23 

0.6402E-02 

0.2076E-02 

91.24 

0.2705E-02 

0.6663E-03 

TT~ 

45.88 

0.1549E-01 

0.5533E-02 

56.02 

0.1019E-01 

0.9948E-02 

61.08 

0.3031E-02 

0.2884E-02 

66.12 

0.7426E-03 

0.9214E-02 

71.17 

0.3778E-03 

0.9359E-03 

76.20 

0.1918E-02 

0.6873E-02 

81.22 

0.6889E-03 

0.1911E-02 

86.23 

0.5036E-02 

0.8409E-02 

91.24 

0.2893E-02 

0.1435E-02 

8.57  MeV  tabulated  data 


^cm 

(deg) 

7r+ 

30.62 

dg 

du/ 

(mb/sr) 

error 

(mb/sr) 

0.3191E-01 

0.6524E-02 

40.80 

0.1193E-01 

0.3340E-02 

45.88 

0.9605E-02 

0.3557E-02 

50.95 

0.3792E-02 

0.1925E-02 

56.02 

0.5158E-02 

0.4972E-02 

61.08 

0.4305E-02 

0.1433E-02 

66.12 

0.3619E-04 

0.1494E-02 

71.17 

0.1262E-02 

0.8460E-03 

76.20 

0.3482E-02 

0.2329E-02 

81.22 

0.7267E-03 

0.7398E-03 

86.23 

0.7550E-03 

0.1442E-02 

91.24 

0.1211E-02 

0.5738E-03 

7r“ 

30.62 

0.1099E-01 

0.8858E-02 

50.95 

0.6287E-02 

0.2530E-02 

66.12 

0.2891E-02 

0.8993E-02 

71.17 

0.3819E-03 

0.9699E-03 

81.22 

0.1884E-02 

0.2042E-02 

86.23 

0.2521E-02 

0.8043E-02 

9.15  MeV  tabulated  data 


^cm 

d<T 

du’ 

error 

(deg) 

-r+ 

(mb/sr) 

(mb/sr) 

TT 

25.53 

0.6742E-01 

0.1252E-01 

30.62 

0.1444E+00 

0.9508E-02 

35.71 

0.1063E+00 

0.7717E-02 

40.80 

0.1137E+00 

0.5771E-02 

45.88 

0.9170E-01 

0.6521E-02 

50.95 

0.6917E-01 

0.3594E-02 

56.02 

0.5291E-01 

0.9790E-02 

61.08 

0.3022E-01 

0.2549E-02 

66.12 

0.2573E-01 

0.3552E-02 

71.17 

0.1807E-01 

0.1788E-02 

76.20 

0.1819E-01 

0.3945E-02 

81.22 

0.2179E-01 

0.1938E-02 

86.23 

0.1862E-01 

0.3165E-02 

91.24 

0.2152E-01 

0.1343E-02 

TT" 

25.53 

0.3199E-01 

0.1664E-01 

30.62 

0.6029E-01 

0.1060E-01 

35.71 

0.8121E-01 

0.1106E-01 

40.80 

0.7659E-01 

0.9921E-02 

45.88 

0.6325E-01 

0.7863E-02 

50.95 

0.4257E-01 

0.3938E-02 

56.02 

0.4230E-01 

0.1331E-01 

61.08 

0.2348E-01 

0.4442E-02 

66.12 

0.5969E-02 

0.1020E-01 

71.17 

0.9399E-02 

0.1648E-02 

76.20 

0.4356E-02 

0.7434E-02 

81.22 

0.9968E-02 

0.3072E-02 

86.23 

0.8556E-02 

0.9289E-02 

91.24 

0.1022E-01 

0.2098E-02 

9.76  MeV  tt*  tabulated  data 


^cm 

error 

(deg) 

(mb/sr) 

(mb/sr) 

TT 

25.53 

0.4064E+00 

0.2430E-01 

30.62 

0.4738E+00 

0.1687E-01 

35.71 

0.4403E+00 

0.1440E-01 

40.80 

0.3422E+00 

0.1044E-01 

45.88 

0.2235E+00 

0.9941E-02 

50.95 

0.1194E+00 

0.4954E-02 

56.02 

0.4080E-01 

0.1027E-01 

61.08 

0.2189E-01 

0.2330E-02 

66.12 

0.1121E-01 

0.2736E-02 

71.17 

0.1712E-01 

0.2133E-02 

76.20 

0.1307E-01 

0.4123E-02 

81.22 

0.1394E-01 

0.1986E-02 

86.23 

0.9110E-02 

0.2502E-02 

91.24 

0.4793E-02 

0.9312E-03 

7r“ 

25.53 

0.5363E+00 

0.3368E-01 

30.62 

0.6783E+00 

0.3011E-01 

35.71 

0.6440E+00 

0.2592E-01 

40.80 

0.4419E-j-00 

0.1960E-01 

45.88 

0.2838E+00 

0.1487E-01 

50.95 

0.1253E+00 

0.5703E-02 

56.02 

0.1086E-01 

0.1022E-01 

61.08 

0.7092E-02 

0.2928E-02 

66.12 

0.2551E-02 

0.1015E-01 

71.17 

0.2116E-01 

0.2321E-02 

76.20 

0.1279E-01 

0.1050E-01 

81.22 

0.1408E-01 

0.4204E-02 

91.24 

0.3701E-02 

0.1528E-02 

9.9  MeV  TT*  tabulated  data 


0cm 

error 

(deg) 

(mb/sr) 

(mb/sr) 

25.53 

0.5195E-01 

0.1892E-01 

30.62 

0.9598E-01 

0.1311E-01 

35.71 

0.4271E-01 

0.9889E-02 

40.80 

0.5256E-01 

0.7517E-02 

45.88 

0.1077E-01 

0.6900E-02 

50.95 

0.9647E-02 

0.3407E-02 

56.02 

0.7046E-02 

0.7473E-02 

61.08 

0.3785E-02 

0.1880E-02 

66.12 

0.1998E-02 

0.2351E-02 

71.17 

0.5099E-02 

0.1710E-02 

76.20 

0.8074E-02 

0.4085E-02 

81.22 

0.4411E-02 

0.1537E-02 

86.23 

0.2994E-02 

0.2309E-02 

91.24 

0.1701E-02 

0.7968E-03 

ir~ 

25.53 

0.2324E+00 

0.2904E-01 

30.62 

0.2242E+00 

0.2483E-01 

35.71 

0.1761E+00 

0.2054E-01 

40.80 

0.1733E+00 

0.1500E-01 

45.88 

0.9126E-01 

0.1245E-01 

50.95 

0.3430E-01 

0.4741E-02 

56.02 

0.1501E-01 

0.1154E-01 

61.08 

0.5336E-02 

0.3036E-02 

66.12 

0.1059E-01 

0.1051E-01 

71.17 

0.7749E-02 

0.2058E-02 

76.20 

0.4113E-02 

0.9650E-02 

81.22 

0.1147E-01 

0.3767E-02 

86.23 

0.1130E-01 

0.9852E-02 

91.24 

0.2563E-02 

0.1484E-02 

10.7  MeV  TT*  tabulated  data 


^cm 

dg 

du/ 

error 

(deg) 

TT'*’ 

25.53 

(mb/sr) 

(mb/sr) 

0.3764E-01 

0.1182E-01 

30.62 

0.2800E-01 

0.7427E-02 

35.71 

0.2281E-01 

0.6013E-02 

40.80 

0.1537E-01 

0.3918E-02 

45.88 

0.2078E-01 

0.4591E-02 

50.95 

0.3128E-01 

0.2933E-02 

56.02 

0.3037E-01 

0.8607E-02 

61.08 

0.3481E-01 

0.2771E-02 

66.12 

0.2906E-01 

0.3983E-02 

71.17 

0.2765E-01 

0.2174E-02 

76.20 

0.3120E-01 

0.4999E-02 

81.22 

0.2011E-01 

0.1945E-02 

86.23 

0.2329E-01 

0.3441E-02 

91.24 

0.1502E-01 

0.1247E-02 

7r“ 

56.02 

0.9440E-02 

0.1073E-01 

61.08 

0.5718E-02 

0.3244E-02 

71.17 

0.3396E-02 

0.1443E-02 

76.20 

0.5127E-03 

0.6930E-02 

81.22 

0.9953E-03 

0.2374E-02 

91.24 

0.3532E-02 

0.1605E-02 

11.3  MeV  TT^  tabulated  data 


^cm 

(deg) 

TT"*" 

35.71 

d<T 

dw 

(mb/sr) 

error 

(mb/sr) 

0.2940E-01 

0.6215E-02 

40.80 

0.9606E-02 

0.3822E-02 

45.88 

0.1223E-01 

0.4406E-02 

50.95 

0.1973E-01 

0.2715E-02 

56.02 

0.1255E-01 

0.6637E-02 

61.08 

0.1310E-01 

0.1828E-02 

66.12 

0.1121E-01 

0.2678E-02 

71.17 

0.6735E-02 

0.1395E-02 

76.20 

0.5845E-02 

0.3143E-02 

81.22 

0.7287E-02 

0.1387E-02 

86.23 

0.7111E-02 

0.2380E-02 

91.24 

0.5348E-02 

0.8875E-03 

7r“ 

40.80 

0.1219E-01 

0.7012E-02 

45.88 

0.3430E-02 

0.5421E-02 

50.95 

0.6317E-02 

0.2997E-02 

56.02 

0.1515E-01 

0.1134E-01 

61.08 

0.1425E-01 

0.3950E-02 

66.12 

0.6439E-02 

0.1087E-01 

71.17 

0.8947E-02 

0.1774E-02 

76.20 

0.4929E-02 

0.7826E-02 

81.22 

0.2613E-02 

0.2540E-02 

91.24 

0.3539E-02 

0.1705E-02 

11.9  MeV  tabulated  data 


^cm 

dff 

du/ 

error 

(deg) 

7r+ 

(mb/sr) 

(mb/sr) 

25.53 

0.1633E+00 

0.1576E-01 

30.62 

0.2151E+00 

0.1147E-01 

35.71 

0.2529E+00 

0.1053E-01 

40.80 

0.1934E+00 

0.7224E-02 

45.88 

0.1290E+00 

0.7901E-02 

50.95 

0.8269E-01 

0.4172E-02 

56.02 

0.5747E-01 

0.1077E-01 

61.08 

0.2254E-01 

0.2559E-02 

66.12 

0.2251E-01 

0.3807E-02 

71.17 

0.1355E-01 

0.1879E-02 

76.20 

0.1380E-01 

0.3d22E-02 

81.22 

0.1256E-01 

0.1714E-02 

86.23 

0.9648E-02 

0.2667E-02 

91.24 

0.1421E-01 

0.1293E-02 

Tt~ 

25.53 

0.3271E+00 

0.2671E-01 

30.62 

0.3598E+00 

0.1811E-01 

35.71 

0.3197E+00 

0.1732E-01 

40.80 

0.2693E+00 

0.1506E-01 

45.88 

0.1970E4-00 

0.1243E-01 

50.95 

0.1078E+00 

0.5814E-02 

56.02 

0.3021E-01 

0.1359E-01 

61.08 

0.3738E-01 

0.5675E-02 

66.12 

0.7687E-02 

0.1176E-01 

71.17 

0.2909E-01 

0.2667E-02 

76.20 

0.1211E-01 

0.9207E-02 

81.22 

0.2232E-01 

0.4400E-02 

86.23 

0.1094E-01 

0.1077E-01 

91.24 

0.1480E-01 

0.2725E-02 

12.5  MeV  TT*  tabulated  data 


Bern 

dg 

error 

(deg) 

tt'*' 

25.53 

(mb/sr) 

(mb/sr) 

0.3588E-01 

0.1242E-01 

30.62 

0.3869E-01 

0.8048E-02 

35.71 

0.5302E-01 

0.6626E-02 

40.80 

0.6132E-01 

0.5230E-02 

45.88 

0.5973E-01 

0.6249E-02 

50.95 

0.6155E-01 

0.3913E-02 

56.02 

0.5233E-01 

0.1076E-01 

61.08 

0.5670E-01 

0.3607E-02 

66.12 

0.5405E-01 

0.5333E-02 

71.17 

0.4594E-01 

0.2940E-02 

76.20 

0.4273E-01 

0.5819E-02 

81.22 

0.3569E-01 

0.2616E-02 

86.23 

0.3557E-01 

0.4200E-02 

91.24 

0.2470E-01 

0.1650E-02 

ir~ 

25.53 

0.1007E+00 

0.2046E-01 

30.62 

0.5242E-01 

0.1098E-01 

35.71 

0.5337E-01 

0.1022E-01 

40.80 

0.4877E-01 

0.8908E-02 

45.88 

0.3199E-01 

0.7357E-02 

50.95 

0.2590E-01 

0.3985E-02 

56.02 

0.2274E-01 

0.1303E-01 

61.08 

0.1416E-01 

0.4396E-02 

66.12 

0.1076E-01 

0.1256E-01 

71.17 

0.4891E-02 

0.1695E-02 

76.20 

0.5812E-02 

0.8844E-02 

81.22 

0.7793E-02 

0.3402E-02 

86.23 

0.1511E-02 

0.9502E-02 

91.24 

0.8373E-02 

0.2290E-02 

12.9  MeV  TT^  tabulated  data 


0cm 

dg 

error 

(deg) 

tt'*' 

30.62 

(mb/sr) 

(mb/sr) 

0.2084E-01 

0.7645E-02 

35.71 

0.1172E-01 

0.7000E-02 

40.80 

0.8090E-02 

0.4489E-02 

45.88 

0.1095E-01 

0.4651E-02 

50.95 

0.8866E-02 

0.2845E-02 

56.02 

0.1980E-01 

0.8385E-02 

61.08 

0.1318E-01 

0.2334E-02 

66.12 

0.1089E-01 

0.3072E-02 

71.17 

0.6324E-02 

0.1576E-02 

76.20 

0.9429E-02 

0.3854E-02 

81.22 

0.3479E-02 

0.1350E-02 

86.23 

0.4098E-02 

0.2468E-02 

91.24 

0.4728E-02 

0.1036E-02 

TT" 

30.62 

0.2092E-01 

0.1023E-01 

35.71 

0.1171E-01 

0.9439E-02 

40.80 

0.5439E-02 

0.7570E-02 

45.88 

0.1253E-01 

0.6010E-02 

50.95 

0.7326E-02 

0.3306E-02 

61.08 

0.9862E-02 

0.4065E-02 

66.12 

0.2583E-02 

0.1151E-01 

71.17 

0.4274E-02 

0.1646E-02 

81.22 

0.8672E-03 

0.2570E-02 

86.23 

0.1664E-02 

0.9965E-02 

91.24 

0.7119E-02 

0.2187E-02 

13.1  MeV  TT^  tabulated  data 


^cm 

da 

error 

(deg) 

tt"*" 

25.53 

(mb/sr) 

(mb/sr) 

0.1858E-01 

0.1185E-01 

30.62 

0.1801E-01 

0.7361E-02 

35.71 

0.6733E-01 

0.8166E-02 

40.80 

0.434VE-01 

0.4970E-02 

45.88 

0.3190E-01 

0.5350E-02 

50.95 

0.1790E-01 

0.2852E-02 

56.02 

0.3087E-03 

0.6636E-02 

61.08 

0.6833E-02 

0.1896E-02 

66.12 

0.2724E-02 

0.2487E-02 

71.17 

0.3152E-02 

0.1306E-02 

81.22 

0.8364E-03 

0.1128E-02 

86.23 

0.3668E-02 

0.2236E-02 

91.24 

0.5685E-02 

0.9567E-03 

7r~ 

25.53 

0.7043E-01 

0.1807E-01 

30.62 

0.4881E-01 

0.1212E-01 

35.71 

0.3377E-01 

0.1064E-01 

40.80 

0.2185E-01 

0.79G2E-02 

45.88 

0.1973E-01 

0.6425E-02 

50.95 

0.1173E-01 

0.3602E-02 

56.02 

0.5593E-02 

0.1213E-01 

61.08 

0.1582E-02 

0.3661E-02 

81.22 

0.4741E-02 

0.2997E-02 

91.24 

0.1723E-02 

0.1893E-02 

14.1  MeV  TT*  tabulated  data 


Bern 

du/ 

error 

(deg) 

25.53 

(mb/sr) 

(mb/sr) 

0.4084E-01 

0.1316E-01 

30.62 

0.6909E-01 

0.9171E-02 

35.71 

0.7331E-01 

0.8340E-02 

40.80 

0.8658E-01 

0.5916E-02 

45.88 

0.8425E-01 

0.7077E-02 

50.95 

0.7711E-01 

0.4221E-02 

56.02 

0.5524E-01 

0.1107E-01 

61.08 

0.5504E-01 

0.3639E-02 

66.12 

0.4467E-01 

0.5039E-02 

71.17 

0.3519E-01 

0.2585E-02 

76.20 

0.2841E-01 

0.5099E-02 

81.22 

0.2718E-01 

0.2364E-02 

86.23 

0.1765E-01 

0.3209E-02 

91.24 

0.2361E-01 

0.1582E-02 

n~ 

25.53 

0.6445E-01 

0.1931E-01 

30.62 

0.1037E+00 

0.1331E-01 

35.71 

0.1321E+00 

0.1288E-01 

40.80 

0.1458E+00 

0.1230E-01 

45.88 

0.1178E+00 

0.1079E-01 

50.95 

0.9559E-01 

0.5804E-02 

56.02 

0.8679E-01 

0.1831E-01 

61.08 

0.5351E-01 

0.6564E-02 

66.12 

0.1663E-01 

0.1431E-01 

71.17 

0.3417E-01 

0.2927E-02 

76.20 

0.1363E-01 

0.1034E-01 

81.22 

0.1738E-01 

0.4405E-02 

86.23 

0.4476E-02 

0.1022E-01 

91.24 

0.1835E-01 

0.3222E-02 

14.4  MeV  7r^  tabulated  data 


^cm 

(deg) 

25.53 

dg 

dw 

(mb/sr) 

error 

(mb/sr) 

0.2939E-01 

0.1249E-01 

30.62 

0.1244E-01 

0.8572E-02 

35.72 

0.2747E-01 

0.7561E-02 

40.80 

0.5308E-01 

0.5845E-02 

45.88 

0.3676E-01 

0.6815E-02 

50.95 

0.2623E-01 

0.3563E-02 

56.02 

0.1463E-01 

0.9010E-02 

61.08 

0.1294E-01 

0.2763E-02 

66.13 

0.3839E-02 

0.3800E-02 

71.17 

0.5695E-02 

0.1697E-02 

76.20 

0.6798E-02 

0.4003E-02 

81.22 

0.5678E-03 

0.1648E-02 

86.24 

0.2898E-02 

0.2682E-02 

91.24 

0.1185E-02 

0.1149E-02 

7r“ 

25.53 

0.4952E-01 

0.1939E-01 

30.62 

0.6839E-02 

0.1155E-01 

35.72 

0.5513E-01 

0.1132E-01 

40.80 

0.4473E-01 

O.lOSOE-01 

45.88 

0.2628E-01 

0.8977E-02 

50.95 

0.3719E-01 

0.4745E-02 

56.02 

0.2352E-01 

0.1G27E-01 

61.08 

0.1105E-01 

0.509GE-02 

66.13 

0.2389E-02 

0.1309E-01 

71.17 

0.3994E-02 

0.2154E-02 

76.20 

0.7401E-02 

0.91G2E-02 

81.22 

0.6277E-02 

0.374GE-02 

86.24 

0.2003E-02 

0.1017E-01 

91.24 

0.1753E-02 

0.2219E-02 

14.6  MeV  TT^  tabulated  data 


^c.„ 

du^ 

error 

(deg) 

TT'*' 

25.53 

(mb/sr) 

(mb/sr) 

0.5434E-01 

0.1295E-01 

30.62 

0.8241E-01 

0.9756E-02 

35.72 

0.7098E-01 

0.8349E-02 

40.80 

0.6560E-01 

0.6468E-02 

45.88 

0.4894E-01 

0.7412E-02 

50.95 

0.3666E-01 

0.3960E-02 

56.02 

0.3590E-01 

0.1055E-01 

61.08 

0.8769E-03 

0.2821E-02 

66.13 

0.1083E-01 

0.37G5E-02 

71.17 

0.5471E-02 

0.2109E-02 

81.22 

0.1017E-01 

0.1940E-02 

86.24 

0.1357E-01 

0.3G42E-02 

91.24 

0.9282E-02 

0.1320E-02 

7r“ 

25.53 

0.5860E-01 

0.19S4E-01 

30.62 

0.9355E-01 

0.137SE-01 

35.72 

0.8263E-01 

0.1227E-01 

40.80 

0.4306E-01 

0.1190E-01 

45.88 

0.6243E-01 

0.1137E-01 

50.95 

0.2929E-01 

0.525GE-02 

56.02 

0.6928E-02 

0.13S()E-01 

61.08 

0.1512E-01 

0.G5.3GE-02 

71.17 

0.1041E-01 

0.2G1SE-02 

81.22 

0.1542E-01 

0.4219E-02 

86.24 

0.1495E-01 

0.1175E-01 

91.24 

0.8580E-02 

0.2572E-02 

15.0  MeV  TT*  tabulated  data 


0cm 

dff 

Alii 

error 

(deg) 

tt'*' 

25.53 

(mb/sr) 

(mb/sr) 

0.9475E-01 

0.1381E-01 

30.62 

0.8344E-01 

0.9147E-02 

35.72 

0.9847E-01 

0.9293E-02 

40.80 

0.1180E+00 

0.7039E-02 

45.88 

0.1203E+00 

0.8811E-02 

50.95 

0.1277E+00 

0.5278E-02 

56.02 

0.9146E-01 

0.1341E-01 

61.08 

0.4757E-01 

0.4397E-02 

66.13 

0.3342E-01 

0.4751E-02 

71.17 

0.1719E-01 

0.2771E-02 

76.20 

0.8628E-02 

0.3998E-02 

81.22 

0.1779E-01 

0.2314E-02 

86.24 

0.2924E-01 

0.4339E-02 

91.24 

0.2706E-01 

0.1862E-02 

7r“ 

25.53 

0.1169E+00 

0.2236E-01 

30.62 

0.1214E+00 

0.1426E-01 

35.72 

0.1377E+00 

0.1401E-01 

40.80 

0.1279E+00 

0.1440E-01 

45.88 

0.1102E+00 

0.1384E-01 

50.95 

0.1242E+00 

0.6740E-02 

56.02 

0.7370E-01 

0.1766E-01 

61.08 

0.8236E-01 

0.1051E-01 

66.13 

0.1674E-01 

0.1390E-01 

71.17 

0.1739E-01 

0.3153E-02 

76.20 

0.3594E-02 

0.8877E-02 

81.22 

0.1733E-01 

0.4458E-02 

86.24 

0.1116E-01 

0.1125E-01 

91.24 

0.2578F-01 

0.3537E-02 

16.5  MeV  TT^  tabulated  data 


^cm 

d<r 

error 

(deg) 

7r+ 

25.53 

(mb/sr) 

(mb/sr) 

0.1117E+00 

0.2550E-01 

30.62 

0.1334E+00 

0.1774E-01 

35.72 

0.8587E-01 

0.1168E-01 

40.80 

0.1184E+00 

0.1140E-01 

45.88 

0.6140E-01 

0.1182E-01 

50.95 

0.4800E-01 

0.6617E-02 

61.08 

0.2757E-01 

0.3989E-02 

66.13 

0.2262E-01 

0.6752E-02 

71.17 

0.2301E-01 

0.3095E-02 

76.20 

0.8762E-02 

0.7493E-02 

81.22 

0.1245E-01 

0.3445E-02 

86.24 

0.2536E-01 

0.6522E-02 

91.24 

0.1625E-01 

0.2699E-02 

TT" 

25.53 

0.1635E+00 

0.4002E-01 

30.62 

0.1068E+00 

0.2417E-01 

35.72 

0.9617E-01 

0.1699E-01 

40.80 

0.4135E-01 

0.1673E-01 

45.88 

0.1955E-01 

0.1325E-01 

50.95 

0.1708E-01 

0.7389E-02 

61.08 

0.4459E-01 

0.8039E-02 

71.17 

0.1716E-01 

0.3108E-02 

86.24 

0.2611E-02 

0.1776E-01 

16.9  MeV  TT*  tabulated  data 


^cm 

d<T 

dti; 

error 

(deg) 

TT"^ 

25.53 

(mb/sr) 

(mb/sr) 

0.1364E+00 

0.1857E-01 

30.62 

0.2051E+00 

0.1288E-01 

35.72 

0.1846E+00 

0.1157E-01 

40.80 

0.1508E+00 

0.8535E-02 

45.88 

0.9650E-01 

0.9440E-02 

50.95 

0.5769E-01 

0.5251E-02 

56.02 

0.3554E-01 

0.1303E-01 

61.08 

0.6247E-02 

0.3503E-02 

66.13 

0.1841E-02 

0.4592E-02 

71.17 

0.1914E-03 

0.2557E-02 

76.20 

0.1629E-01 

0.6547E-02 

81.22 

0.6389E-02 

0.2699E-02 

86.24 

0.4206E-02 

0.3299E-02 

91.24 

0.9425E-02 

0.1949E-02 

7r“ 

25.53 

0.2458E+00 

0.2763E-01 

30.62 

0.2355E+00 

0.1849E-01 

35.72 

0.2068E+00 

0.1885E-01 

40.80 

0.1286E+00 

0.1488E-01 

45.88 

0.8546E-01 

0.1211E-01 

50.95 

0.2607E-01 

0.5883E-02 

56.02 

0.5142E-02 

0.1589E-01 

61.08 

0.1800E-02 

0.6539E-02 

66.13 

0.5427E-02 

0.1650E-01 

71.17 

0.1449E-02 

0.2918E-02 

76.20 

0.1315E-02 

0.1194E-01 

81.22 

0.5061E-02 

0.4488E-02 

86.24 

0.7475E-02 

0.1306E-01 

91.24 

0.7396E-02 

0.3597E-02 

17.2  MeV  7r^  tabulated  data 


^cm 

dw 

error 

(deg) 

(mb/sr) 

(mb/sr) 

TT 

25.53 

0.1007E+00 

0.1783E-01 

30.62 

0.7009E-01 

0.1088E-01 

35.72 

0.6831E-01 

0.1033E-01 

40.80 

0.7131E-01 

0.7755E-02 

45.88 

0.6541E-01 

0.1105E-01 

50.95 

0.4674E-01 

0.6145E-02 

56.02 

0.7418E-01 

0.1534E-01 

61.08 

0.7432E-01 

0.5819E-02 

66.13 

0.7319E-01 

0.6859E-02 

71.17 

0.6217E-01 

0.4667E-02 

76.20 

0.6797E-01 

0.9079E-02 

81.22 

0.6787E-01 

0.4268E-02 

86.24 

0.6025E-01 

0.6001E-02 

91.24 

0.3641E-01 

0.2569E-02 

7C~ 

25.53 

0.1067E+00 

0.2434E-01 

30.62 

0.9298E-01 

0.1641E-01 

35.72 

0.7468E-01 

0.1725E-01 

40.80 

0.1056E+00 

0.1430E-01 

45.88 

0.6970E-01 

0.1447E-01 

50.95 

0.7932E-01 

0.7054E-02 

56.02 

0.5613E-01 

0.2070E-01 

61.08 

0.1094E+00 

0.1007E-01 

66.13 

0.9792E-01 

0.2324E-01 

71.17 

0.1046E+00 

0.6042E-02 

76.20 

0.6389E-01 

0.1761E-01 

81.22 

0.8812E-01 

0.8888E-02 

86.24 

0.5544E-01 

0.1705E-01 

91.24 

0.4593E-01 

0.5227E-02 

17.6  MeV  TT*  tabulated  data 


^cni 

du^ 

error 

(deg) 

(mb/sr) 

(mb/sr) 

25.53 

0.1608E+00 

0.2043E-01 

30.62 

0.2003E+00 

0.1531E-01 

35.72 

0.1493E+00 

0.1364E-01 

40.80 

0.1542E+00 

0.1003E-01 

50.95 

0.4770E-01 

0.6207E-02 

56.02 

0.3834E-01 

0.1504E-01 

61.08 

0.1072E-01 

0.5948E-02 

66.13 

0.1839E-01 

0.5820E-02 

71.17 

0.1471E-01 

0.4748E-02 

76.20 

0.2635E-01 

0.7688E-02 

81.22 

0.2848E-01 

0.3615E-02 

86.24 

0.2791E-01 

0.5775E-02 

91.24 

0.3022E-01 

0.2559E-02 

7r“ 

25.53 

0.2554E4-00 

0.3393E-01 

30.62 

0.2301E+00 

0.2184E-01 

35.72 

0.2503E+00 

0.2207E-01 

40.80 

0.1698E+00 

0.1895E-01 

50.95 

0.5982E-01 

0.7685E-02 

56.02 

0.5743E-01 

0.2337E-01 

66.13 

0.1112E-01 

0.2124E-01 

71.17 

0.2066E-01 

0.5536E-02 

76.20 

0.3863E-01 

0.1810E-01 

81.22 

0.2646E-01 

0.7206E-02 

86.24 

0.2466E-01 

0.1768E-01 

91.24 

0.2562E-01 

0.4672E-02 

18.3  MeV  tabulated  data 


^cm 

du; 

error 

(deg) 

;r+ 

25.53 

(mb/sr) 

(mb/sr) 

0.1091E+00 

0.1649E-01 

30.62 

O.llOlE+00 

0.1184E-01 

35.72 

0.6236E-01 

0.1262E-01 

40.80 

0.8940E-01 

0.7697E-02 

45.88 

0.6781E-01 

0.8294E-02 

50.95 

0.3015E-01 

0.4394E-02 

61.08 

0.6283E-02 

0.3343E-02 

66.13 

0.1643E-01 

0.4679E-02 

71.17 

0.8873E-02 

0.2540E-02 

76.20 

0.1786E-01 

0.5777E-02 

81.22 

0.1312E-01 

0.2579E-02 

86.24 

0.6570E-02 

0.3747E-02 

91.24 

0.6486E-02 

0.1583E-02 

TT~ 

25.53 

0.1452E+00 

0.2589E-01 

30.62 

0.1006E+00 

0.1592E-01 

35.72 

0.9299E-01 

0.1666E-01 

40.80 

0.4642E-01 

0.1643E-01 

45.88 

0.6658E-01 

0.1221E-01 

50.95 

0.2151E-01 

0.5624E-02 

61.08 

0.1327E-01 

0.6142E-02 

66.13 

0.4440E-02 

0.1582E-01 

71.17 

0.1670E-01 

0.2875E-02 

76.20 

0.1134E-01 

0.1222E-01 

81.22 

0.1803E-01 

0.5285E-02 

86.24 

0.1816E-01 

0.1473E-01 

91.24 

0.6354E-02 

0.4034E-02 

18.7  MeV  tabulated  data 


Ocn 

du; 

error 

(<leg) 

(mb/sr) 

(mb/sr) 

TT 

25.53 

0.4253E-01 

0.1321E-01 

30.62 

0.6870E-01 

0.9777E-02 

35.72 

0.2240E-01 

0.9525E-02 

40.80 

0.3670E-01 

0.6315E-02 

45.88 

0.4168E-01 

0.7520E-02 

50.95 

0.2426E-01 

0.4241E-02 

56.02 

0.1304E-01 

0.1165E-01 

61.08 

0.6683E-02 

0.2981E-02 

66.13 

0.8795E-02 

0.4455E-02 

71.17 

0.1273E-01 

0.2487E-02 

76.20 

0.1245E-01 

0.5065E-02 

81.22 

0.5254E-02 

0.2326E-02 

86.24 

0.1280E-01 

0.4317E-02 

91.24 

0.9816E-02 

0.1922E-02 

TT" 

25.53 

0.1107E+00 

0.2433E-01 

30.62 

0.1036E+00 

0.1518E-01 

35.72 

0.6292E-01 

0.1332E-01 

40.80 

0.6864E-01 

0.1445E-01 

45.88 

0.1659E-01 

0.9594E-02 

50.95 

0.1053E-01 

0.4781E-02 

66.13 

0.2982E-02 

0.1339E-01 

71.17 

0.7664E-02 

0.2405E-02 

76.20 

0.1160E-01 

0.1146E-01 

81.22 

0.7798E-02 

0.5418E-02 

86.24 

0.3451E-02 

0.1220E-01 

91.24 

0.4988E-02 

0.2726E-02 

18.9  MeV  TT^  tabulated  data 


0cm 

du; 

error 

(deg) 

(mb/sr) 

(mb/sr) 

TT 

25.53 

0.1692E-01 

0.1300E-01 

30.62 

0.2407E-01 

0.9072E-02 

35.72 

0.3340E-01 

0.9364E-02 

40.80 

0.1855E-01 

0.6251E-02 

45.88 

0.3824E-01 

0.7344E-02 

50.95 

0.1872E-01 

0.4235E-02 

56.02 

0.1035E-01 

0.1210E-01 

61.08 

0.1381E-01 

0.3167E-02 

66.13 

0.2452E-01 

0.5195E-02 

71.17 

0.9216E-02 

0.2460E-02 

76.20 

0.1056E-01 

0.5536E-02 

81.22 

0.1801E-01 

0.2700E-02 

86.24 

0.1164E.01 

0.4334B-02 

91.24 

0.5471E-02 

0.1875E-02 

7r“ 

25.53 

0.5522E-01 

0.2280E-01 

30.62 

0.3946E-01 

0.1507E-01 

35.72 

0.6037E-01 

0.1313E-01 

40.80 

0.1638E-01 

0.1358E-01 

45.88 

0.2042E-01 

0.9659E-02 

50.95 

0.1129E-01 

0.4849E-02 

61.08 

0.9398E-02 

0.5738E-02 

66.13 

0.2444E-02 

0.1334E-01 

71.17 

0.7298E-02 

0.2384E-02 

81.22 

0.1164E-01 

0.5669E-02 

86.24 

0.1344E-01 

0.1331E-01 

91.24 

0.7229E-02 

0.3175E-02 
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